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Introductory Remarks. 

JL HE former Essay* contained a general method for reducing all the most important 
problems of dynamics to the study of one characteristic function, one central or ra- 
dical relation. It was remarked at the close of that Essay, that many eliminations 
required by this method in its first conception, might be avoided by a general trans- 
formation, introducing the time explicitly into a part S of the whole characteristic 
function V ; and it is now proposed to fix the attention chiefly on this part S, and to 
call it the Principal Function. The properties of this part or function S, which were 
noticed briefly in the former Essay, are now more fully set forth ; and especially its 
uses in questions of perturbation, in which it dispenses with many laborious and cir- 
cuitous processes, and enables us to express accurately the disturbed configuration of 
a system by the rules of undisturbed motion, if only the initial components of veloci- 
ties be changed in a suitable manner. Another manner of extending rigorously to 
disturbed motion the rules of undisturbed, by the gradual variation of elements, in 
number double the number of the coordinates or other marks of position of the 
system, which was first invented by Lagrange, and was afterwards improved by 
Poisson, is considered in this Second Essay under a form perhaps a little more ge- 
neral ; and the general method of calculation which has already been applied to 
other analogous questions in optics and in dynamics by the author of the present 
Essay, is now applied to the integration of the equations which determine these ele- 
ments. This general method is founded chiefly on a combination of the principles of 
variations with those of partial differentials, and may furnish, when it shall be ma- 
tured by the labours of other analysts, a separate branch of algebra, which may be 
called perhaps the Calculus of Principal Functions ; because, in all the chief applica- 
tions of algebra to physics, and in a very extensive class of purely mathematical 
questions, it reduces the determination of many mutually connected functions to the 
search and study of one principal or central relation. When applied to the integration 
of the equations of varying elements, it suggests, as is now shown, the consideration 

* Philosophical Transactions for the year 1834, Second Part. 
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of a certain Function of Elements, which may be variously chosen, and may either 
be rigorously determined, or at least approached to, with an indefinite accuracy, 
by a corollary of the general method. And to illustrate all these new general 
processes, but especially those which are connected with problems of perturbation, 
they are applied in this Essay to a very simple example, suggested by the motions of 
projectiles, the parabolic path being treated as the undisturbed. As a more important 
example, the problem of determining the motions of a ternary or multiple system, 
with any laws of attraction or repulsion, and with one predominant mass, which was 
touched upon in the former Essay, is here resumed in a new way, by forming and inte- 
grating the differential equations of a new set of varying elements, entirely distinct 
in theory (though little differing in practice) from the elements conceived by La- 
grange, and having this advantage, that the differentials of all the new elements for 
loth the disturbed and disturbing masses may be expressed by the coefficients of one 
disturbing function. 

Transformations of the Differential Equations of Motion of an Attracting or Repelling 

System. 

1 . It is well known to mathematicians, that the differential equations of motion of 
any system of free points, attracting or repelling one another according to any func- 
tions of their distances, and not disturbed by any foreign force, may be comprised in 
the following formula : 

2. w(3"&ar+y&y + *"&*) =&U: (1.) 

the sign of summation 2 extending to all the points of the system ; m being, for any 
one such point, the constant called its mass, and wyz being its rectangular coordi- 
nates ; while x"y"z" are the accelerations, or second differential coefficients taken 
with respect to the time, and ix, ty, I z are any arbitrary infinitesimal variations of 
those coordinates, and U is a certain force-function, introduced into dynamics by La- 
grange, and involving the masses and mutual distances of the several points of the 
system. If the number of those points be ?i, the formula (1.) may be decomposed into 
3 n ordinary differential equations of the second order, between the coordinates and 
the time, 

and to integrate these differential equations of motion of an attracting or repelling 
system, or some transformations of these, is the chief and perhaps ultimately the only 
problem of mathematical dynamics. 

2. To facilitate and generalize the solution of this problem, it is useful to express 
previously the 3 n rectangular coordinates xy z as functions of 3n other and more 
general marks of position *ii*i 2 ... n 3n ; and then the differential equations of motion 
take this more general form, discovered by Lagrange, 
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d 8T 8T _ 8U 
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in which 

T=±2\m(at 2 +y' 2 + z' 2 ) 

For, from the equations (2.) or (1.), 
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in which 



and 
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T being here considered as a function of the 6 n quantities of the forms tf and n, ob- 
tained by introducing into its definition (4.), the values 

Ix 
"Ji 



i i S x i i 8 x i i / 

* =^8T + J? 28^ + ---+"s«8, 



"Ja 



-, &c. 



'Sn 



(8.) 



A different proof of this important transformation (3.) is given in the Meeanique 
Analytique. 

3. The function T being homogeneous of the second dimension with respect to the 
quantities q', must satisfy the condition 

2T = 2-«W-> • • • (»•) 

and since the variation of the same function T may evidently be expressed as follows, 

ST = 2(|fw + ^), (10.) 

we see that this variation may be expressed in this other way, 

>T=2(,/>$-£>„) (11.) 

If then we put, for abridgement, 
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ST 
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and consider T (as we may) as a function of the following form, 

T = F (w l9 Wg, • • • ™ 3n , *li, *l2> • • • %r), ■ • • • 



we see that 



and 



SF 



= »i, 



8F 



8F 
8 1i 



ST 



= n 



8F 
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(13.) 
(14.) 



ST 



and therefore that the general equation (3.) may receive this new transformation, 

d*r. __ S(U-F) 
~di ~ *\ 



(15.) 



(16.) 



If then we introduce, for abridgement, the following expression H, 



H = F - U = F(w x , w 2} . . . vr 3n , n x , n 2 , . • .fi 3t ) -U(j? 15 ?? 2 , ,. . n Sn ), . (17.) 

we are conducted to this new manner of presenting the differential equations of 
motion of a system of n points, attracting or repelling one another : 
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dt 
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(A.) 



In this view, the problem of mathematical dynamics, for a system of n points, is to 
integrate a system (A.) of 6 n ordinaiy differential equations of the first order, be- 
tween the 6m variables n i **. and the time t ; and the solution of the problem must 
consist in assigning these 6m variables as functions of the time, and of their own 
initial values, which we may call e. p.. And all these 6n functions, or Qn relations 
to determine them, may be expressed, with perfect generality and rigour, by the 
method of the former Essay, or by the following simplified process. 

Integration of the Equations of Motion, by means of one Principal Function, 
4. If we take the variation of the definite integral 

S =X( 2 -*^- H ) rf ' (18 * } 

without varying t or dt, we find, by the Calculus of Variations, 

i§ = r t i$.dt, (19.) 

«/ o 



in which 



6 ATT 3 
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and therefore 

*S' = 2(^H-*H^), (2L) 

that is, by the equations of motion (A.), 

SS' = 2(^ + ^) = ^2.^ ; (22.) 

the variation of the integral S is therefore 

§ S = 2 (sr In - p i e), (23.) 

(p and e being still initial values,) and it decomposes itself into the following 6 n ex- 
pressions, when S is considered as a function of the 6 n quantities ^. <?., (involving also 
the time,) 

ss ss 

w i = H" ; Pl ~ ~ ^l ; 






r > ...•....".... (B.) 



SS __ _8S 



''3ra "S n 

which are evidently forms for the sought integrals of the 6 n differential equations of 
motion (A.), containing only one unknown function S. The difficulty of mathema- 
tical dynamics is therefore reduced to the search and study of this one function S, 
which may for that reason be called the Principal Function of motion of a system. 
This function S was introduced in the first Essay under the form 

the symbols T and U having in this form their recent meanings ; and it is worth 
observing, that when S is expressed by this definite integral, the conditions for its 
variation vanishing (if the final and initial coordinates and the time be given) are 
precisely the differential equations of motion (3.), under the forms assigned by La- 
grange. The variation of this definite integral S has therefore the double property, 
of giving the differentia! equations of motion for any transformed coordinates when 
the extreme positions are regarded as fixed, and of giving the integrals of those dif- 
ferential equations when the extreme positions are treated as varying. 

5. Although the function S seems to deserve the name here given it of Principal 
Function, as serving to express, in what appears the simplest way, the integrals of the 
equations of motion, and the differential equations themselves ; yet the same analy- 
sis conducts to other functions, which also may be used to express the integrals of 
the same equations. Thus, if we put 

«=X(- 2 -'W+n)^, <*j 

and take the variation of this integral Q without varying t or dt, we find, by a simi- 
lar process, 

>Q = 2(?U-e^); (25.) 

o 2 
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so that if we consider Q as a function of the 6 n quantities sr. p. and of the time, we 
shall have 6 n expressions 

"i = + 8^.' e i = ~ 8^7> ( 26 -) 

which are other forms for the integrals of the equations of motion (A.), involving the 
function Q instead of S. We might also employ the integral 

8 H „ . « r» 



V 



=X 2 -<*W^ = 2 ./>^> (27.) 

which was called the Characteristic Function in the former Essay, and of which, when 
considered as a function of the 6 n + 1 quantities j?. e. H, the variation is 

&V = 2 (w &«-/»& e)+*&H (28.) 

And all these functions S, Q, V, are connected in such a way, that the forms and 
properties of any one may be deduced from those of any other. 

Investigation of a Pair of Partial Differential Equations of the first Order, which the 

Principal Function must satisfy. 

6. In forming the variation (23.), or the partial differential coefficients (B.), of the 

Principal Function S, the variation of the time was omitted ; but it is easy to calcu- 

s s 
late the coefficient -jj corresponding to this variation, since the evident equation 

dS 8S , _.'*S dri . v 

77 = 77 + 2 77 77 • ■ • (29-> 

gives, by (20.), and by (A.), (B.), 

I7 = S'-2. w ^r= -H. (30.) 

It is evident also that this coefficient, or the quantity — H, is constant, so as not 
to alter during the motion of the system ; because the differential equations of mo- 
tion (A.) give 

dH _ /SH d n 8Hrf<*\ 

dt — Z \ir, dt + S« dt) ~ U> ' • • * • ' ' • " • ( 31 -) 

If, therefore, we attend to the equation (17-)) an d observe that the function F is neces- 
sarily rational and integer and homogeneous of the second dimension with respect to 
the quantities wr., we shall perceive that the principal function S must satisfy the two 
following equations between its partial differential coefficients of the first order, 
which offer the chief means of discovering its form : 

8S /*S 8S 8S \ TT .1 

17 + F VH ' &V ■ • ' 8^' "i» '» • • • %«) = U (*i» >l» • • • * 3 J, | 

8S /8S 8S 8S \ TT . 

77 + F (t7;> 87"' • ' • 87~' e i> e 2, • • • * 3n ) = U (e l5 e 2 , . . . e 3n ). 



>(C) 
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Reciprocally, if the form of S be known, the forms of these equations (C.) can be 
deduced from it, by elimination of the quantities e or y between the expressions of its 
partial differential coefficients ; and thus we can return from the principal function S 
to the functions F and U, and consequently to the expression H, and the equations 
of motion (A.). 

Analogous remarks apply to the functions Q and V, which must satisfy the partial 
differential equations, 
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and 






sv sv 
sv sv 

S^ ' ' ' * Sc ' 



• • • «3n) = H + U (j? l5 »» 2 , . . . ? 8 J, | 



e i> e 2 , . . . e 3B ) = H + U (e u e 2 , . . . e 3n ). 



(33.) 



General Method of improving an approximate Expression for the Principal Function 

in any Problem of Dynamics. 

7. If we separate the principal function S into any two parts, 

S 1 + S 2 = S, (34.) 

and substitute their sum for S in the first equation (C), the function F, from its 
rational and integer and homogeneous form and dimension, may be expressed in this 
new way, 

„/8S SS 

F ls 



>■ (35.) 



because 



and 



,sv • • 8, 3 ;^i» • -W - * W • • *v 3n ' ni ' ' * %n > 

^/SS, SS, \ „/8S 2 SS 2 \ 

+ r«5)iJ + ... + r(E)& 

^^)=^H)-F'fti:), o«.) 



2 . f (t) 



SS, 



and since, by (A.) and (B.), 
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Tit I 8 S\ _,, , . 8H dri. . ' 

F k ^) = p to = nr, = ir ............. (38.) 

we easily transform the first equation (C.) to the following, 

§=-^+u(» l ,..o-F(if;,.|^,, 1 ,., 3n )+p(^.-^» I .-'s.)> < D > 

which gives rigorously 

^-jC{-^+ u (*-%J-'C^-t^.*.-*.)}"] 



(E.) 



■>2n 



supposing only that the two parts S 1} S 2 , like the whole principal function S, are 
chosen so as to vanish with the time. 

This general and rigorous transformation offers a general method of improving an 
approximate expression for the principal function S, in any problem of dynamics- 
For if the part Sj be such an approximate expression, then the remaining part S 2 will 
be small ; and the homogeneous function F involving the squares and products of the 
coefficients of this small part, in the second definite integral (E.), will be in general 
also small, and of a higher order of smallness ; we may therefore in general neglect 
this second definite integral, in passing to a second approximation, and may in general 
improve a first approximate expression S x by adding to it the following correction, 

AS 1 =X{-^ + U ^--0-F(|f J ..|A,, 1 ,.., 3ra )}^ ; (R) 

in calculating which definite integral we may employ the following approximate forms 
for the integrals of the equations of motion, 

8S, SSj 8S, , , 

Pi = -8^2 =-^-,..^3* =-8^' • • • ( 39 -> 

* * on 

expressing first, by these, the variables % as functions of the time and of the 6 n con- 
stants ej p b and then eliminating, after the integration, the 3 n quantities p h by the same 
approximate forms. And when an improved expression, or second approximate value 
S 1 + A S l5 for the principal function S, has been thus obtained, it maybe substituted 
in like manner for the first approximate value S i} so as to obtain a still closer ap- 
proximation, and the process may be repeated indefinitely. 

An analogous process applies to the indefinite improvement of a first approximate 
expression for the function Q or V. 

Rigorous Theory of Perturbations, founded on the Properties of the Disturbing Part 

of the whole Principal Function. 

8. If we separate the expression H (17.) into any two parts of the same kind, 

H X + H 2 =H, (40.) 
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In which 

H x = F x fo, w 2 , . . . ar 3n , ??i, % • • "s n) — U l (*i> "25 • • >?s«)> (41.) 

and 

H 2 = F 2 (*!, « a , . . . w 3 „, 9i, * 2 > • • O — U 2 fe *2> • • *s»)» ..... (42.) 

the functions F x F 2 U x U 2 being such that 

F 1 +F 2 = F,U 1 + U 2 = U; (43.) 

the differential equations of motion (A.) will take this form, 

15 - iMi a. IMs — { - *Jk_*Jk /m 

dt ~ Sot. + Sw.' dt — H S«. J l"-J 

and if the part H 2 and its coefficients be small, they will not differ much from these 
other differential equations, 

so that the rigorous integrals of the latter system will be approximate integrals of 
the former. Whenever then, by a proper choice of the predominant term H x , a 
system of 6 n equations such as (H.) has been formed and rigorously integrated, 
giving expressions for the 6 n variables jj £ m i as functions of the time t, and of their 
own initial values e { p which may be thus denoted : 

ii= : 'Pi(t>e 1 ,e 2 ,..e 3n ,p 1 ,p 2 ,..p3n), (44.) 

and 

Vf i = '4'i(*,e 1 ,e 2 ,..e 3n ,p 1 ,p 2 ,..p Sn ); (45.) 

the simpler motion thus defined by the rigorous integrals of (H.) may be called the 
undisturbed motion of the proposed system of n points, and the more complex motion 
expressed by the rigorous integrals of (G.) may be called by contrast the disturbed 
motion of that system ; and to pass from the one to the other, may be called a Pro- 
blem of Perturbation. 

9. To accomplish this passage, let us observe that the differential equations of un- 
disturbed motion (H.), being of the same form as the original equations (A.), may 
have their integrals similarly expressed, that is, as follows : 

SS t 8S t ,, v 

V *-J7* Pi = -17.* (L - ) 

Sj being here the principal function of undisturbed motion, or the definite integral 

s i=X( 2 --S- H 0^ • < 46 -) 

considered as a function of the time and of the quantities % e,- In like manner if we 
represent by S x + S 2 the whole principal function of disturbed motion, the rigorous 
integrals of (G.) may be expressed by (B.), as follows : 

"i- 17 + 17> Pi~~ 17- 17. ^ 
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Comparing the forms (44.) with the second set of equations (I.) for the integrals of 
undisturbed motion, we find that the following relations between the functions <p t Sj 
must be rigorously and identically true : 

% = «Pi^e lje2 ,.. e3n ,-g-,-g-,..-g— J; ....... (47.) 

>■ *■ oil 

and therefore, by (K.), that the integrals of disturbed motion may be put under the 

following forms, 

/ . ■ 8S 2 . 8S 2 ■ 8S 2 \ 

V i = p i \t,e 1 ,e 2 ,..e 3n ,p l + j-,p 2 + j-,..p 3n + sj-). . . . . (L.) 

1 2 3n 

We may therefore calculate rigorously the disturbed variables m by the rules of un- 
disturbed motion (44.), if without altering the time /, or the initial values e { of those 
variables, which determine the initial configuration, we alter (in general) the initial 
velocities and directions, by adding to the elements p t the following perturbational 
terms, 

A ^ = lf> A P2 = |fV- Aj» 3ll = ^: ......... (M.) 

a remarkable result, which includes the whole theory of perturbation. We might 

deduce from it the differential coefficients j/j 5 or the connected quantities zj { , which 

determine the disturbed directions and velocities of motion at any time t ; but a 

similar reasoning gives at once the general expression, 

8 S 2 , / ■ 8 S 2 8 S 2 8 So \ ._ T 

™i ~ J7 + 4i {*, *i> *v • • e sn»Pi + TZ> Pz + T7> • • • P*» + YTJ* • ( N -> 

implying, that after altering the initial velocities and directions or the elements p. as 
before, by the perturbational terms (M.),we may then employ the rules of undisturbed 
motion (45.) to calculate the velocities and directions at the time t, or the varying 
quantities &., if we finally apply to these quantities thus calculated the following new 

corrections for perturbation : 

a 8 S 2 8 S 2 . 8 So 

Aw i = V Aw 2 = h»-- a ^ = ^ ......... (O.) 

Approximate expressions deduced from the foregoing rigorous Theory. 

10. The foregoing theory gives indeed rigorous expressions for the perturbations, 
in passing from the simpler motion (H.) or (I.) to the more complex motion (G.) or 
(K.) : but it may seem that these expressions are of little use, because they involve an 
unknown disturbing function S 2 , (namely, the perturbational part of the whole princi- 
pal function S,) and also unknown or disturbed coordinates or marks of position ?/.. 
However, it was lately shown that whenever a first approximate form for the princi- 
pal function S, such as here the principal function S x of undisturbed motion, has been 
found, the correction S 2 can in general be assigned, with an indefinitely increasing 
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accuracy ; and since the perturbations (M.) and (O.) involve the disturbed coordi- 
nates n. only as they enter into the coefficients of this small disturbing function S 2 , it 
is evidently permitted to substitute for these coordinates, at first, their undisturbed 
values, and then to correct the results by substituting more accurate expressions. 

1 1 . The function S x of undisturbed motion must satisfy rigorously two partial dif- 
ferential equations of the form (C), namely, 
8S, 



j 

it 



+ F l UV • • • 0~^> "I'-- *S«) = U l 0l> • • • V3n)> 

8 S x -j-, /8 S t 8 S t \ . 



(P-) 



and therefore, by (D.),the disturbing function S 2 must satisfy rigorously the following 
other condition : 

~^f =U 2 ( ni , . . n sn )- F 2 ^, . . g— -, n lt . . , 3fl j + F^, . . g— , , l5 . . « 3n ), (Q.) 

and may, on account of the homogeneity and dimension of F, be approximately ex- 
pressed as follows : 

S 2 = Jl { U 2 ("1' ' • ^3 J - F 2 (g^ • • §^S "I, • • 1 3 J } dt > • • • ( R 

or thus, by (I.), 

S *-X { U 2^1' * '«sJ ~ F 2 K> • • "s.»«l» • -"sj} «**> .... (S.) 

that is, by (42.), 

S 2 =~f o l Yl 2 dt (T.) 

In this expression, H 2 is given immediately as a function of the varying quantities 
q. to-., but it may be considered in the same order of approximation as a known func- 
tion of their initial values e. p. and of the time t, obtained by substituting for n to. 
their undisturbed values (44.) (45.) as functions of those quantities; its variation 
may therefore be expressed in either of the two following ways : 

SH 2 = 2(^^+^Sto), (48.) 

or 

U% = ?(^le + ^t>p)+^U (49.) 

Adopting the latter view, and effecting the integration (T.) with respect to the 
time, by treating the elements e. p. as constant, we are afterwards to substitute for 
the quantities p. their undisturbed expressions (39.) or (I.), and then we find for the 
variation of the disturbing function S 2 the expression 

*S 2 =-n 2 Zt + 2(-hef;^d t + l^.£^dt), . (50.) 

mdcccxxxv. p 
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which enables us to transform the perturbational terms (M.) (O.) into the following 
approximate forms : 

m„ ,. . „ PS, /*'8H 3 



A / H8H *,7,_l_V PS ' /* 



$ p 



dt, 



and 



Aw. = 2 






(U.) 



(V.) 



8e8» i «/o 8p 

containing only functions and quantities which may be regarded as given, by the 
theory of undisturbed motion. 

12. In the same order of approximation, if the variation of the expression (44.) for 
an undisturbed coordinate t). be thus denoted, 

the perturbation of that coordinate may be expressed as follows : 



Sp 



that is, by (U.), 

8», ajh 



Ap; 



8«. 



(W.) 



-* 8 H, 



8 ^ S 2 Si 



k / $p at 

*8H. 



3re 



"t- V8 Pl 8^ T 8p 2 te t 8e 2 T ' ' ' + 8 p 3n ie^e^JJo fp~ x rf ' 



+ 



}> • (52.) 






8 a Sj 8,, i_J!A. 



+ 



Besides, the identical equation (47.) gives 

H __ H _f_Sj_ H _8^Si_ , 
Se k ■ip l 9e t te l + Sp a Se k ie i ~ f ' 



+ ...+ 



H 8 2 S! 



8» Q 8e 3 Q 

•« 3 re 3 re 

8 1* 8 a S t 



)/: 






dt. 



" ' 8 ^ 3 re H»3re 

the expression (52.) may therefore be thus abridged, 



(53.) 



8«. /»i 

A*. = -— J / 

8m. 



8H„ 



rf£- 



*P3nJ<> 

8 », 



*8H S 
8> 



rf£ 



3n 



^l^ 8 ^l 8 's*«A> »*„„ ' J 



(X.) 



and shows that instead of the rigorous perturbational terms (M.) we may approxi- 
mately employ the following, 

**= -£'-£**** (YO 

l 

in order to calculate the disturbed configuration at any time t by the rules of undis- 
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turbed motion, provided that besides thus altering the initial velocities and directions 
we alter also the initial configuration, by the formula 

^=jf^ ( Z -) 

It would not be difficult to calculate, in like manner, approximate expressions for the 
disturbed directions and velocities at any time t ; but it is better to resume, in an- 
other way, the rigorous problem of perturbation. 

Other Rigorous Theory of Perturbation, founded on the properties of the disturbing 
part of the constant of living force, and giving formulae for the Variation of Ele- 
ments more analogous to those already known. 

13. Suppose that the theory of undisturbed motion has given the 6n constants 
e. p. or any combinations of these, » v » 2 , . . . » gw , as functions of the 6n variables 

n. m. and of the time t, which may be thus denoted : 

and which give reciprocally expressions for the variables q. ■&. in terms of these ele- 
ments and of the time, analogous to (44.) and (45.), and capable of being denoted 
similarly, 

*i = Pi (*> *1> *2» • • • *6»)> ™i = ^i &> Xi > %2 > ' • * *6f) ; • • • • ( 55< ) 

then, the total differential coefficient of every such element or function x., taken with 

respect to the time, (both as it enters explicitly and implicitly into the expressions 
(54.),) must vanish in the undisturbed motion ; so that, by the differential equations 
of such motion (H.), the following general relation must be rigorously and identically 
true: 

= T i+2( T -*i5!_.8H!) (56 . 

In passing to disturbed motion, if we retain the equation (54.) as a definition of the 
quantity *, that quantity will no longer be constant, but it will continue to satisfy 

the inverse relations (55.), and may be called, by analogy, a varying element of the 
motion ; and its total differential coefficient, taken with respect to the time, may, by 
the identical equation (56.), and by the differential equations of disturbed motion 
(G.), be rigorously expressed as follows : 

^= 2 (!5?L&_i5iH,) (A ,) 

at \«i| ow 6 -as o »j / 

14. This result (A 1 .) contains the whole theory of the gradual variation of the ele- 
ments of disturbed motion of a system ; but it may receive an advantageous trans- 
formation, by the substitution of the expressions (55.) for the variables n. vs. as func- 
tions of the time and of the elements ; since it will thus conduct to a system of 6 n 

p2 



108 PROFESSOR HAMILTON ON A GENERAL METHOD IN DYNAMICS. 

rigorous and ordinary differential equations of the first order between those varying 
elements and the time. Expressing, therefore, the quantity H 2 as a function of these 
latter variables, its variation & H 2 takes this new form, 

SH 2 =2.^* + 8 -#^, (57 .) 

and gives, by comparison with the form (48.), and by (54.), 

8H 2 _ ^ SH 8 h 8H 3 _ ^ 8H 2 8* 
8» _i, b 8i ' $u _Z ' h b ; • • • • • • ( bH -> 

ftp '.j» v> j* 

and thus the general equation (A 1 .) is transformed to the following, 



d* t 8H, S_H 8_H, 

* ■* on 



in which 

(8x. 8 >t 8 k. 8x \ 
17 "8^ ~ "s^ 1^7 : ( C ^ 

so that it only remains to eliminate the variables y\ w from the expressions of these 
latter coefficients. Now it is remarkable that this elimination removes the symbol t 
also, and leaves the coefficients a is expressed as functions of the elements x alone, not 
explicitly involving the time. This general theorem of dynamics, which is, perhaps, 
a little more extensive than the analogous results discovered by Lagrange and by 
Poisson, since it does not limit the disturbing terms in the differential equations of mo- 
tion to depend on the configuration only, may be investigated in the following way. 

15. The sign of summation 2 in (C 1 .), like the same sign in those other analogous 
equations in which it has already occurred without an index in this Essay, refers not 
to the expressed indices, such as here i, s, in the quantity to be summed, but to an 
index which is not expressed, and which may be here called r ; so that if we intro- 
duce for greater clearness this variable index and its limits, the expression (C 1 .) be- 
comes 

3« / 8x. 8x 8x. 8x \ , 

a is = 2 (•'_*'): (59.) 

and its total differential coefficient, taken with respect to the time, may be separated 
into the two following parts, 



+ ( f) i \8« p i/ H iv r dtSi r )> 



> (60.) 



which we shall proceed to calculate separately, and then to add them together. By 
the definition (54.), and the differential equations of disturbed motion (G.), 
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in which, by the identical equation (56.), 



8 3 X £ g 3n /8x 

8* Sot. Sot, (m) 



» /^ 8_Hj H SJHA _ 

i \S « Sot 8ot 8 m / ' 
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(62.) 



we have therefore 






d 

dt 8 or 

3 B / 8 X.. 



8H 2 

J5f 8W 



S 2 x. »n 8x 

2 w -tio . i 



8ot 8i 



+ 



■H, 



hi. it 'u r 



__» i!iLi_ 

H, 8 ^„&ot 



), • (63.) 



and 

8«/JX. ,7 X 
y lid s 

l V 8m dt 8ot 



may be found from this, by merely changing itos: so that 

8 x Ox £ 8 x; 



r) = 



5m d£ S-st 
8 2 x. 8 x ; 8 2 x. 



-8 x. 8 2 x 



3n,3nj-/8x, 8' x, _ 8 x. ^ x g H / x . a* », x & 0< x. y Hj 

(r,«)i,i I \o % ov u ov r 8« r 8 WK 8W 8m u +\8« r !,,8 Wf 8,, r 8^ Sot,./ 8ot m 

+ V8, r 8 Wfi 8^/ S^Sot." 1 " Vs.,,.8^ ir> r t\) ** n *«J ' 



\ (64.) 



and similarly 

3n / 8 x ,7 8 x. 
I s a i 



8x. 



8 x. 8 2 x 



8x 8 2 x. 



8H, 



(65.) 



(r)l\t*r r dt X 0K r dt OYlJ ' 

3 n,sn c , 8x, _8^_ __ 8x. l*x x g H2 /8*. J^ _ ^x_ s 6-x. X 
(r,.)i,i IVSotS^S^ 8ot S^Vs^ "^VSot 8^8,, Sot^S^ 

+ VSot 8« m 8ot 8« m ; Sot^ "*" W, Sot m Sot Sot/ 8„ m o % J ' 

Adding, therefore, the two last expressions, and making the reductions which pre- 
sent themselves, we find, by (60.), 



d „ 3» / .(« 

t-: a. =2 A 



(«)8H, , D (»)8H 
* r -t> "» — : 

OH ' ,• ,, OT 



:)■ 



(D».) 



in which 
A = 



(«) 
B = 



3n / 8x 8" J X. 

2 ( 



+ 



8x. SV 



8 x 8 3 x. 



(r)i \8.) r 8OT M SOT r 8i) p 8w tt 8«r p 8w r 8ot w 8^ 8ot ? , 8OT M 8t, r . 



;). 



3n / 8x 8 a x. 
(r)l \SOT p S\Sl) ( . 



t / 

2 ■(: 



8x. 8 s x 



8ot 8 yi 8ti 



+ 



8x. 



8x 8 s x. 



8ti 8 « 8 vi 

'r 'II r 



K 8 '. 8 



i«r / 



(66.) 



and since this general form (D 1 .) for j- f a i>s contains no term independent of the dis- 

8 H 8 H 
turbing quantities — -, -y- 3 , it is easy to infer from it the important consequence 

already mentioned, namely, that the coefficients a i>s , in the differentials (B 1 .) of the 
elements, may be expressed as functions of those elements alone, not explicitly in- 
volving the time. 
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It is evident also, that these coefficients a it s have the property 

%* = *-%*, (670 

and 

%t = 0; . (68.) 

the term proportional to y^ 2 disappears therefore from the expression (B 1 .) for-—; 

and the term 

S H 2 8 H 2 . 8 H 2 dx; 

a, „ . -*— - in 



8x, * "'.« " 8x, '" 8X; rf< ' 



"-s 



destroys the term 

8H a 8H g . 8H 2 dx, 

8x s ,a «.** 8^ ln 8x, dt ' 

when these terms are added together ; we have, therefore, 

*-'-£lB-°. < E '-) 

or 

tf* ~ 8* ' (*'•) 

that is, in taking the first total differential coefficient of the disturbing expression H 2 
with respect to the time, the elements may be treated as constant. 

Simplification of the differential equations which determine these gradually varying 
elements, in any problem of Perturbation ; and Integration of the simplified equations 
by means of certain Functions of Elements. 

16. The most natural choice of these elements is that which makes them corre- 
spond, in undisturbed motion, to the initial quantities e i p t . These quantities, by the 
differential equations (H.), may be expressed in undisturbed motion as follows, 

ei=: " i ~X^7 dt} * = •* +X M? ** ; (69 - } 

and if we suppose them found, by elimination, under the forms 

e i ~n i + ®i (t, *!, y 2 , . . . y 3n , ar lt * 2 , . . . *f 3n ), 1 

Pi = *| + % (t, til, %, • • • n Zn , *i, * 2 » • • • *Sn)» J 

it is easy to see that the following equations must be rigorously and identically true, 
for all values of n- t w 

= ** (0, ti v y 2 , ... n$ n , «i, *a, . . . * 8B ), 1 

= T, (0, jj 15 ?? 2 , . . . ?? 3n , Wj, «r 2 , . . . *,„).■ J 

When, therefore, in passing to disturbed motion, we establish the equations of defi- 
nition, 
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»l— Vi+ G*«(Ml>*2» •••*8»>*1»*2» ■••• or S »)» 1 , „ 

f (72.) 

introducing 6 n varying elements » t \, of which the set \ would have been represented 
in our recent notation as follows : 

we see that all the partial differential coefficients of the forms _J, _J, _J, *, vanish 

oij r 8«r r 8ij r Sw r 

when £ = 0, except the following : 

~ = 1, ^ — 1 ; (74.) 

and, therefore, that when t is made = 0, in the coefficients a i} „ (59.), all those coeffi- 
cients vanish, except the following : 

a r,3n f r = ** a Sn + r,r = ~ 1 (75.) 

But it has been proved that these coefficients a if „ when expressed as functions of 
the elements, do not contain the time explicitly ; and the supposition t = introduces 
no relation between those 6 n elements *,. \, which still remain independent : the co- 
efficients a i>s , therefore, could not acquire the values 1, 0, —1, by the supposition 
t = 0, unless they had those values constantly, and independently of that supposition. 
The differential equations of the forms (B 1 .), may therefore be expressed, for the pre- 
sent system of varying elements, in the following simpler way : 

~dt ~ 8a. ; ~di- 8*. ' V* * 

and an easy verification of these expressions is offered by the formula (E 1 .), which 
takes now this form, 

2(^ + ^w)=° (H.) 

17- The initial values of the varying elements *,- \ L are evidently e i p i) by the defi- 
nitions (72.), and by the identical equations (71.) ; the problem of integrating rigo- 
rously the equations of disturbed motion (G.), between the variables r\ { zs t and the 
time, or of determining these variables as functions of the time and of their own 
initial values e { p t , is therefore rigorously transformed into the problem of integrating 
the equations (G 1 .), or of determining the 6 n elements x, h Xi as functions of the time 
and of the same initial values. The chief advantage of this transformation is, that if 
the perturbations be small, the new variables (namely, the elements,) alter but little : 
and that, since the new differential equations are of the same form as the old, they 
may be integrated by a similar method. Considering, therefore, the definite integral 
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as a function of the time and of the 6n quantities » v * 2 , . . . z 3n , e x , e 2 , . . . e 3n , and 
observing that its variation, taken with respect to the latter quantities, may be shown 
by a process similar to that of the fourth number of this Essay to be 

&E = 2(U*-/>&e), (I 1 .) 

we find that the rigorous integrals of the differential equations (G 1 .) may be ex- 
pressed in the following manner : 

o E 5 E , ' . 

\ = T7.> Pi ~-T7.> • ( K -) 

■i i 

in which there enters only one unknown function of elements E, to the search and 
study of which single function the problem of perturbation is reduced by this new 
method. 
We might also have put 

C=^(-2.*^ + H 2 )^, .(77.) 

and have considered this definite integral C as a function of the time and of the 6 n 
quantities X { p { ; and then we should have found the following other forms for the in- 
tegrals of the differential equations of varying elements, 

I S C 8 C ,T 1 X 

** = + IX.' e i = - 1¥. • • • • ( L •) 

i •*» 

And each of these functions of elements, C and E, must satisfy a certain partial differ- 
ential equation, analogous to the first equation of each pair mentioned in the sixth 
number of this Essay, and deduced on similar principles. 

18. Thus, it is evident, by the form of the function E, and by the equations (K 1 .), 
(G 1 .), and (76.), that the partial differential coefficient of this function, taken with 
respect to the time, is 

8J_rfE 8E d* _ 

and therefore that if we separate this function E into any two parts 

Ej + E^E, (N<.) 

and if, for greater clearness, we put the expression H 2 under the form 

H 2 = H 2 (t, * 1} * 2 , . . . » 3n) Xj, x 2 , . . . x 3n ), (O 1 .) 

we shall have rigorously the partial differential equation 

O-^.SE , S IE, m+i^X. fp M 

V — It + It "r **2 [*»*!»" -*3n»s x + 8->---8jl +8x Q )• V •> 

\ I ' 3 n 3 ii/ 

which gives, approximately, by (G 1 .) and (K 1 .), when the part E 2 is small, and when 
we neglect the squares and products of its partial differential coefficients, 

d Jn + Hi . H (t x , 1A 8EA 



= 



E 2 
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Hence, in the same order of approximation, if the part E l5 like the whole function E, 
be chosen so as to vanish with the time, we shall have 

= -X{^+h 2 (m..-^-©} • • ■ • • (*'■> 

and thus a first approximate expression E x can be successively and indefinitely cor- 
rected. 
Again, by (L 1 .) and (G 1 .), and by the definition (77-), 

8 C dC ^SCrfX __ -e,, v 

U — 1t~ 2 • J\ ~dl ~ H 2 ; • • • • • • ( & ■) 

the function C must therefore satisfy rigorously the partial differential equation, 

L2 _ h ( t — — > > V <rn 

and if we put 

C = C 1 + C 2J (IF.) 

and suppose that the part C 2 is small, then the rigorous equation 

8/ "•" H — n 2 V r ' 8x, T" 8V'" "8\, "T~ Sa ' / 'i'" ,A 3V • " ' * ^ v '> 

1 * 3 re 3 m 

becomes approximately, by (G 1 .) and (L 1 .), 

dt ~ It ^ "2 V>tK 1 >---%\ 3n > K l>--- A 3nJ> \ yy •> 

and gives by integration, 

c 2 =x{- s # + H 2 ( ( ,|a > ...|i ) x 1; ...x„,)}^, . . . . pp.) 

the parts C x and C 2 being supposed to vanish separately when t = 0, like the whole 
function of elements C. 

And to obtain such a first approximation, Ej or C l5 to either of these two functions 
of elements E, C, we may change, in the definitions (76.) (77-), the varying elements 
*, X, to their initial values e, p, and then eliminate one set of these initial values by 
the corresponding set of the following approximate equations, deduced from the for- 
mulae (G 1 .) : 



«* = *«+/ T?<**; " ' " • <YM 



and 



h P . 



= Pi-XTt dt (ZK) 



It is easy also to see that these two functions of elements C and E are connected 
with each other, and with the disturbing function S 2 , so that the form of any one 
may be deduced from that of any other, when the function S L of undisturbed motion 
is known. 
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Analogous formulae for the motion of a Single Point. 

19. Our general method in dynamics, though intended chiefly for the study of 
attracting and repelling systems, is not confined to such, but may be used in all 
questions to which the law of living forces applies. And all the analysis of this 
Essay, but especially the theory of perturbations, may usefully be illustrated by the 
following analogous reasonings and results respecting the motion of a single point. 

Imagine then such a point, having for its three rectangular coordinates x y z, and 
moving in an orbit determined by three ordinary differential equations of the second 
order of forms analogous to the equations (2.), namely, 

f7 _ SU „_8U „_8U 
x ~ %x> y ~ §y> Z ~ 8z ' ' ' \'°'' 

U being any given function of the coordinates not expressly involving the time : and 

let us establish the following definition, analogous to (4.), 

T = ~(*'2+y2 + *'V • (79.) 

x' y' z' being the first, and x" y" z" being the second differential coefficients of the 
coordinates, considered as functions of the time t. If we express, for greater gene- 
rality or facility, the rectangular coordinates x y zas functions of three other marks of 
position q x n 2 j? 3 , T will become a homogeneous function of the second dimension of 
their first differential coefficients n\ ii% V3 taken with respect to the time ; and if we 
put, for abridgement, 

___ 8T __ ST __ ST 

OT i — Wi' ™ 2 ~ Wf Ws "~ *7? ' (80.) 

T mav be considered also as a function of the form 

T = F (ra-j, otvj, vr 3 , n x , ti 2 , i} 3 ), (81.) 

which will be homogeneous of the second dimension with respect to nt x sr 2 nr 3 . We 
may also put, for abridgement, 

F (ro-j, w 2 , «r 3 , m, n 2 , tf 3 ) — U (fli, n 2 , n 3 ) = H ; . . (82.) 

and then, instead of the three differential equations of the second order (78.), we may 
employ the six following of the first order, analogous to the equations (A.), and ob- 
tained by a similar reasoning, 



dt~-" i "iv l } dt~~ t l^' dt * "*" Sw 3 ' I 
dia l 8 H dvs % 8H diiT 3 8 H 



(83.) 



dt ~ &V dt 8>j 2 > dt S >j 3 " 

20. The rigorous integrals of these six differential equations may be expressed 
under the following forms, analogous to (B.), 

8S 8S 8S ■) 



?V "2-8,8' w 3~ g^ I 



05;,,- " 0)].,- i 

(84.) 



8S _ 8_S 8S 
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in which e x e 2 e 3 p 1 p 2 p 3 are the initial values, or values at the time 0, of m % j? 3 
tETj rs 2 «r 3 ; and S is the definite integral 

S =/"(•> if + «*{f + -»I5l-H)", (85.) 

considered as a function of m *i% % e x e 2 e 3 and t. The quantity H does not change 
in the course of the motion, and the function S must satisfy the following pair of 
partial differential equations of the first order, analogous to the pair (C), 

8S . -/8S 8S SS \ TT , . 

"87 + F W W **' *' * 2 ' " 3 / = U ( " 15 * 2 ' * s) ; 



8S 
8* 



, -n/8S 8S 8S \ TT , N 



This important function S, which may be called the principal function of the motion, 
may hence be rigorously expressed under the following form, obtained by reasonings 
analogous to those of the seventh number of this Essay : 

+y f Vh-h'k-H'H-8^'^^^^; j 

S x being any arbitrary function of the same quantities ^ v\ 2 n ti e 1 e 2 e 3 t, so chosen as 
to vanish with the time. And if this arbitrary function S t be chosen so as to be a 
first approximate value of the principal function S, we may neglect, in a second ap- 
proximation, the second definite integral in (870 • 

21. A first approximation of this kind can be obtained, whenever, by separating 
the expression H, (82.), into a predominant and a smaller part, Hx and H 2 , and by 
neglecting the part H 2 , we have changed the differential equations (83.) to others, 
namely, 



(88.) 



and have succeeded in integrating rigorously these simplified equations, belonging to 
a simpler motion, which may be called the undisturbed motion of the point. For the 
principal function of such undisturbed motion, namely, the definite integral 

S i=X(-i|f +^ + -3^-H 1 )^,. (89.) 

considered as a function of ri x q 2 q 3 e x e 2 e 3 t, will then be an approximate value for the 
original function of disturbed motion S, which original function corresponds to the 
more complex differential equations, 

Q2 



rfij, _8H, ^ S _8H, dvi 3 8H t 
dt " 8w, 5 dt ~~ 8w 2 ' dt " 8«r 3 ' 




1 

r 

j 


dvr x 8 Hj d ot 2 8 H t d «r 3 
dt ~ " Sjjj' dt ~~ ' 8ij 2 ' dt ~~ 


8H! 

8 Is' 
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dt 



— LSl .IS <Hz_ 8H, SHg <H,_ SHi 



rf Wj 8 Hj 



8 H 2 <£ «r 2 
TX/' ~dT 



8H, 

8 1 2 



8 H 2 <i! Wg 

1X/>~d7 : 



8 H, 8 H 2 



(90.) 



'% 



8 % 'J 



The function S x of undisturbed motion must satisfy a pair of partial differential 
equations of the first order, analogous to the pair (86.) ; and the integrals of undis- 
turbed motion may be represented thus, 



'i — 



8S, 
8V 



8Sj 



8X> OT 3 



8_S, 
8%' 



8 S! 8 Sj 8 S[ 

Pi = ~ IX' -^2 — ~ 8X' 7^3 IT" : 



> 



(91.) 



8^' - - „ C2 - - vc 3 

while the integrals of disturbed motion may be expressed with equal rigour under the 
following analogous forms, 

8Sj . 8S 9 SS, . 8S 2 8S, , 8S 2 1 



m i~ 8,, + SV 



8 S t , 8 S 2 _ 8 S t , 8 S 2 

r 2 ~ 8X + Ttf OT 3 - g^j- + g^= 



P\ = 



8S, 
8*i 



8_S 2 
8<?i' #* 



8 S, 8 So 



8e„ 



8<? 2 ' ^3— 8g 3 



8S t __ 8_So | 



\ 



(92.) 



if S 2 denote the rigorous correction of S l5 or the disturbing part of the whole principal 
function S. And by the foregoing general theory of approximation, this disturbing 
part or function S 2 may be approximately represented by the definite integral (T.), 

S 2 = — J H 2 d * > - • * 



(93.) 



in calculating which definite integral the equations (91.) may be employed. 
22. If the integrals of undisturbed motion (91.) have given 
Vi = Pi {t, e lt e 2 , e 3 , p v p 2 , p 3 ), 

*2 = P2 (*» e V e 2> e 3> Pv i>2> Pz)> 

% = P 3 if, e v e 2 , e 3 , Pi,p 2 ,p- 6 ), 



(94.) 



and 



1 



w l = ^1 ('» e l> e 2> e i> PV P2> Ps)> 

<*2 = ^2 (<> e l> e 2> e V PV ?2> Ps)> / ■ 

w 3 = ^3 (*> e V e 2> e 3> PV P& Pi) J 



(95.) 



then the integrals of disturbed motion (92.) maybe rigorously transformed as follows, 

8 S„ . 8 So . 8 S 



- / , 1 8 bo . 8 bo ,8 b,\ 

*1 = Pi \*> e V e 2> e 3>Pl + jX> Pi + J^> Pi + jfj, 

*. ( j. 1 8S 8 8 So . 8S 2 \ 

= ?2 {?> e l> e 2' ^ Pi + f^,P2 + J% P A + jfj, > 

*. ( 1 1 8 So . 5S, ,8 So\ 

= Ps \t, e» e 2 , e A , Pl + j^, p 2 + g -£, p 3 + jfj, 



n 2 



h 



(96.) 
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and 

8 So , , / , . 8S 8 . 8 S 2 ,8 S 2 \ 

w i = g^ + ^i \*> e » e * ^i + J7S v* + 8v #» + r^> 

8 S 2 . , / , , 8 S 9 8 S 2 .8 S 2 \ . 

^2 = ^f + ^ 2 (/> e i5 e 25 <& V\ + f^, Pi + if^ Ps + iy s )> I * ( 97 ") 

8S 2 , , / , 8S 2 , 8S 2 , 8S 2 \ 

w 3 = 8^ + ^3 (', *i> «fc ^Pi + sv ^2 + 8^' ^ + g^ > j 

S 2 being here the rigorous disturbing function. And the perturbations of position, at 
any time t, may be approximately expressed by the following formula, 

1 teiJo 8 i>i § Wo 0^2 s Wo 8 j°s I 

_|i/»»i^_|i/«^a^_Ja.r»L^ f ' (98 ° 

8 iWo 8<? i 8 iWo 8e a 8 iWo 8 * 3 J 

together with two similar formulae for the perturbations of the two other coordinates, 
or marks of position q 2 , fj z . In these formulae, the coordinates and H 2 are supposed 
to be expressed, by the theory of undisturbed motion, as functions of the time t, and 
of the constants e x e 2 e 3 p 1 p 2 p 3 . 

23. Again, if the integrals of undisturbed motion have given, by elimination, ex- 
pressions for these constants, of the forms 

e i = "i + 3»i (t> "u n 2 , %, »d ®2, %), "j 

«2 = "2 + #2 ( f > "l5 "25 %> W l, ^2, W 3 ), >........ (99). 

e i = % + ^3 ('j "l5 "25 "35 W 15 ^25 & 3)> J 

and 

J»l = OT 1 + "*1 (*> "l5 "25 "35 ^15 w 25 ™z), ") 

J»2 = ^2 + "**2 (^5 "l5 "25 "35 ^15 ^25 ^3)5 >••...... (100.) 

I 

^3 = ^3 + ¥3 (*5 "15 "25 "35 W 15 OT 25 OT 3) 5 J 

and if, for disturbed motion, we establish the definitions 

X l = "l + ^1 ('» "lJ "25 "35 ™1, ^25 W 3)5 "j 

*2 = "2 + 3*2 (^5 "15 "25 "35 ^15 W 25 ^3)5 /»........ (101.) 

*3 = "3 + *3 (*5 "15 "25 "35 1*1, ™2, ^3)3 J 

and 

h = "*>! + "*1 ('> "I' "25 "35 m l, ™2, ^3)5 1 

^2 = w 2 + ^2 ( f , "15 "25 "35 ^15 w 25 ^3)5 r" • (102.) 

^3 ~ ^3 + ^3 ('> "15 "25 "35 OT D ^25 ^3) » J 

we shall have, for such disturbed motion, the following rigorous equations, of the 
forms (94.) and (95.), 

"l = Pi {*, H, H, »35 ^15 hf h)> I 

"2 = ?2 ( f > x i> x 2, x 3> h, h, h), > .......... (103.) 

"3 = ?3 (^ x v x 2, x 3, h, h, h), J 
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and 



*1 — ^l {t> *U *2> *3> \) ^2> ^3) 

5 

I 



:1 



ro 2 = ^ 2 (/, * l9 *2, * 3 , X u X 2 , X 3 ), }> (104.) 

W 3 == ^3 (*> K l> x 2> x 3> h> h> h) ' J 

and may call the quantities x x x 2 * 3 X x X 2 X 3 the 6 varying elements of the motion. To 
determine these six varying elements, we may employ the six following rigorous equa- 
tions in ordinary differentials of the first order, in which H 2 is supposed to have been 
expressed by (103.) and (104.) as a function of the elements and of the time : 

dx l _ 8 H 2 dx% _ 8H 2 dxz _ 8 H 2 -^ 

dt ~~ 8a,' dt ~ 8a> dt ~ Sa„' 1 

^ 3 l (105.) 

ih —. _ ilk ^? - __ Lli ^3 _ _ LMs . ! 

dt ~ 8x,'» d* ~~ 8x 2 > -d< — 8x3 ' J 

and the rigorous integrals of these 6 equations may be expressed in the following 

manner, 

8E 8E _8E 

X i — 8x,» ^ — 8x,» X 3 — 8x a ' 

8E 8E »* <* • * * * ' ' (106>) 

ft =: ~ K' ft = ~ W ft = 
the constants e x e 2 e 3 p x p 2 p 3 retaining their recent meanings, and being therefore the 
initial values of the elements x x x 2 x 2 h h h > while the function E, which may be 
called the function of elements, because its form determines the laws of their variations, 
is the definite integral 




E 



considered as depending on x x x 2 « 3 e x e 2 e 3 and/. The integrals of the equations (105.) 

may also be expressed in this other way, 

8C S_C 8C 1 

«i — + j » » *2 — + 8 A > *3 — + g a> J 

>....... (108.) 

8C 8C 8C V ' 

e i - ~ 8 A » e 2 - ~ 8p 2 > e 3 - - 8ft » J 
C being the definite integral 

regarded as a function of ^ X 2 X 3 p x p 2 p 3 and £ : and it is easy to prove that each of 
these two functions of elements, C and E, must satisfy a partial differential equation 
of the first order, which can be previously assigned, and which may assist in disco- 
vering the forms of these two functions, and especially in improving an approximate 
expression for either. All these results for the motion of a single point, are analogous 
to the results already deduced in this Essay, for an attracting or repelling system. 
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Mathematical Example, suggested by the motion of Projectiles. 

24. If the three marks of position m fy % °f tne moving point are the rectangular 
coordinates themselves, and if the function U has the form 

Urr-^-il/^CV + ^ + 'V}* (HO.) 

g, p, v, being constants ; then the expression 

H = iK2 + ^2 + ra . 3 2)+g, 3 + |{|«. 2 K 2 + »? 2 2 )+v2^}, . . (HI.) 

is that which must be substituted in the general forms (83.), in order to form the 6 
differential equations of motion of the first order, namely, 

dni <hs dla „ j 

c • • • 



dt 

dw l 



dn*_ drja 

■i» dt """ r 8» dt 



»3> 






d «*» 



2 dvf z _ 



vj 



(112.) 



ST-"!* *> IT 
These differential equations have for their rigorous integrals the six following, 



7} 

^ = e x cos w- 1 -f- — sin p t, 

r* 
Tf 

n 2 = e 2 cos (h t + — sin p £, 

^ 3 = e$ cos c f + — sin p £ — 4- vers n *, 



> • 



(113.) 



and 



1 



(114.) 



w 1 = j9 x cos (&t — ft* e x sin ju. £, 
cr 2 = p 2 cos pi — (fie 2 sin jia £, 

«r 3 = /» 3 cos » t — I v e 3 + -7) sin c * ; 

e i e 2 e 3 Pi P2.P3 heing still the initial values of ^ n % 9) 3 vr t «r 2 ^ 

Employing these rigorous integral equations to calculate the function S, that is, by 
(85.) and (110.) (111.), the definite integral 

8 =J«(lL±j£±*£ + v)dt, (115.) 

we find 

i( Wl 2+V + ^ 2 ) = i{pi 2 + P2 2 +P3 2 + f* 2 K 2 + e 2 2) + (v e3 + f) 2 jl 

+ |{Pi 2 +P2 2 -i«' 2 (ei 2 + ^2 2 )}cos2 / */-|^(^ 1 + e 2 ^ 2 )sin2p/ >(116.) 



and 



U = 



H- ifa - (»<% + f) 2 }cos2^ - i(vc 3 + f ) ft sin2»f, 
{pi 2 + P2 2 + Pi 2 + t* 2 (ei 2 + *2 2 ) + (>e 3 + f) 2 } 



1 



+ i {?i 2 + P2 2 - I* 2 (<?i 2 -I- ea 2 )} cos2f*# - }p (e^i + e^ft) sin 2 f* f f (117-) 
+ i {pi- (»«3 + f) 9 ]cos2v* - I (v% + ■f)p 3 sin2»*j 
and therefore, 
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s 



_£* 



2 v s 



+ iPi 2 +J»2 2 ~ i» 2 (*i 2 + e 2 2 )} Si2 ~f i - i (^ + e 2 p 2 ) vers2^"j 



+ {^-(- 3 +f)7^^-^3(e3 + f)ver S 2^. 



>(118.) 



In order, however, to express this function S,as supposed by our general method, in 
terms of the final and initial coordinates and of the time, we must employ the analo- 
gous expressions for the constants^ p 2 p z , deduced from the integrals (113.), namely, 
the following : 






ft ijj — fj. e x cos ju. t 
sin fit 

^ >i 2 — /x e 2 cos ft t 
sin ft, t 

g 



"13 + 



("* + f) 



» e 3 + -2- 1 cos v ( 



smvt 



and then we find 

Q _ tl . SI 
2 V 2 " 1 " 2 



(119.) 



tanfU ~^ 2 



(l 3 ~ g 3 ) 2 
tan v£ 



1 



- p fa fj + n 2 e 2 ) tan -^ i> (^ + 7) (e 3 + 7) tan ^. 



(120.) 



(121.) 



This principal function S satisfies the following pair of partial differential equations 
of the first order, of the kind (86.), 

8* + 2 (u^J + \8^/ + V8>j a / j - #*3 aWi+12; 2^3, 

ss , 1 r/8s\ 3 . /8s\ s . /ts\n ■ ^ , 2 . 2N * 2 2 

17 + ¥{W +W + Vs^j j = -^ e 3-^(^ 2 + ^ 2 )-^^ 2 ; 

and if its form had been previously found, by the help of this pair, or in any other 
way, the integrals of the equations of motion might {by our general method) have been 
deduced from it, under the forms, 



zs-y — g — = p fa — e{) cotan p t — (h e x tan -g-, 
sr 2 = s — = (J* fa ~~ e 2) cotan pi — pe 2 tan -3—, 



w 3 = §T = » Ofc ~ «s) cotan »f - (»e 3 + 7) 



and. 



tan 



^ 



(122.) 



joj = — g — = («- (^ — ej cotan ^t -\- ^vi v tan -^- , 
^2= — sT" = l* i(«a ~ ^2) cotan /«/# + ^jj 2 tan -g-, 



(123.) 



B = - 87; = » Ok - e 3 ) cota » » « + ("» 3 + -7) tan ^ : 

the last of these two sets of equations coinciding with the set (119.), or (113.), and 
conducting, when combined with the first set, (122.), to the other former set of inte- 
grals, (114.). 
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25. Suppose now, to illustrate the theory of perturbation, that the constants p, v 
are small, and that, after separating the expression (111.) for H into the two parts, 

H 1 = ^(V + ^2 2 + ^3 2 ) + S'^ (124.) 



and 



H 9 



i 2 - H/* 2 (fl 2 + >?2 2 ) + * 2 *3 2 }> (125.) 

we suppress at first the small part H 2 , and so form, by (88.), these other and simpler 
differential equations of a motion which we shall call undisturbed: 



(126.) 



dt 


— ro"u 


dt 


— 


W 2 5 


dt 


= 


OT 3J 


d ■ar l 
dt 


= 0, 


dts^ 
dt 




o, 


du s 
dt 


— 


~~&' - 



These new equations have for their rigorous integrals, of the forms (94.) and (95.) 
>Ji = e 1 +p 1 t, ti 2 — e 2 +p 2 t, % = e 3 +p 3 t— ^gt 2 , . . 



and 



"l=Pl> ™2—P2> * f 3 = P3 



gt; 



(127.) 
(128.) 



and the principal function S t of the same undisturbed motion is, by (89.), 

s i =Jo \- — i — - e %) dt 



=X (*' + *' +ft ' ~ge 3 - 2gp 3 t+g*t>) dt ►, 

= fi'+ft' + ft- -^e 3 ) * - ^ 3 ' 2 + i^ 2 ^ 
or finally, by (127-), 

J>1 — ST 2 6 * W3 + e 3-> 24 © * ' 



(129.) 



(130.) 



This function satisfies, as it ought, the following pair of partial differential equa- 
tions, 

8S X . 1 r^SA 2 , ASA 2 , /SSA 2 1 



Si! 

8_S 
8/ 



•'+i{Gt)*+Gt)'+(Jf)>-^ 



(131.) 



And if by the help of this pair, or in any other way, the form (130.) of this principal 
function S t had been found, the integral equations (127-) and (128.) might have been 
deduced from it, by our general method, as follows : 



r - . 



*2 = g 



8Sj _ l! — e t 

8 i)i ~" t > 

8 Sj la ~ ■ e i 

8 Sl _ lj 8 - g 3 _ X 



*s-W 3 -~ 



gt> 



(132.) 



MDCCCXXXV. 



R 



122 

and 



PROFESSOR HAMILTON ON A GENERAL METHOD IN DYNAMICS. 



fli 



SS, 



8 Sj % — g g 

#2 — ~ FT — 7 



P3--F7: 



% 



(133.) 



the latter of these two sets coinciding with (127), and the former set conducting to 
(128.). 

26. Returning now from this simpler motion to the more complex motion first 
mentioned, and denoting by S 2 the disturbing part or function which must be added 
to S x in order to make up the whole principal function S of that more complex mo- 
tion ; we have, by applying our general method, the following rigorous expression 
for this disturbing function, 

in which we may, approximately, neglect the second definite integral, and calculate 
the first by the help of the equations of undisturbed motion. In this manner we find, 
approximately, by (125.) (127-), 

-H 2 =-^ {(e l + Pl t)»+(e 2 +p 2 ty }-£(et+p 3 t-igfi)*, (135.) 

and therefore, by integration, 



S 2 — - ¥ { ^ ( e l 2 + ^) + " 2 *3 2 > ' - "« {^ ( e lPl + e 2P2) + »* ^Ps) t 2 
-T{^(Pl 2 +P2 2 )+" 2 (P, 2 -g^}f+^^gp^-^gU^ 



1 



(136.) 



or, by (133.), 



s 2 = - %- (V + e \ *i + e \ + n£ + e%n 2 -\- <? 2 2 ) ] 

-^{e + e3*3 + e3 2 +Tg^ + %) t2 + Tog 2t *} : J 



(137.) 



the error being of the fourth order, with respect to the small quantities p, v. And 
neglecting this small error, we can deduce, by our general method, approximate forms 
for the integrals of the equations of disturbed motion, from the corrected function 
S x -f S 2 , as follows : 

S Sn . 8 S 2 r ll -e l p» t ( ,1 \ 

8 Si . 8 S 2 >) 2 - <? 2 \>?t ( ,1 \ 

™* = Ti + ji : =-r 1 --TK»2 + %*%)> >• (138.) 

8S! , 8 8, ij 3 -e 3 1 x »«* / ,1 , 1 A 
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and 



- __L§,_8_s 2 _ !!L - il , £LL . 

Pl~ le x ie l ~ t + 3 l e ' + 



■Vi 



__ 8'S 2 ijg — g a ,(^t( I J_ \ 

^2 — - g^ 2 - Y^ — ~~1 r 3 V e 2 -r 2 V* 

or, in the same order of approximation, 

"i = e i + Pi* - "a ^ 2 * 2 ( e i + jft *)' 

^2 = e 2+?2 # ~ a" ^ 2 ^ 2 ( e 2 + 1"P2V' "* 

i .„ i „.„/ .1 • __* ^ 



(139.) 



and 



OT i = ^i - P 2 t (e x + g- />! *), 

^—P2 — ^ 2 t (e 2 + -%p 2 t), ; 



(140.) 



(141.) 



Accordingly, if we develope the rigorous integrals of disturbed motion, (113.) and 
(114.), as far as the squares (inclusive) of the small quantities ^ and v, we are con- 
ducted to these approximate integrals ; and if we develope the rigorous expression 
(120.) for the principal function of such motion, to the same degree of accuracy, we 
obtain the sum of the two expressions (130.) and (137.)- 

27- To illustrate still further, in the present example, our general method of suc- 
cessive approximation, let S 3 denote the small unknown correction of the approximate 
expression (137-)? so that we shall now have, rigorously, for the present disturbed 
motion, 

S = S L + S 2 + S 3 , (142.) 

Sj and S 2 being here determined rigorously by (130.) and (13/.). Then, substituting 
Sj -4- S 2 for Sj in the general transformation (87.)> we find, rigorously, in the present 
question, 






(143.) 



and if we neglect only terms of the eighth and higher dimensions with respect to 
"the small quantities p, v, we may confine ourselves to the first of these two definite 
integrals, and may employ, in calculating it, the approximate expressions (140.) for 



r2 
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the coordinates of disturbed motion. In this manner we obtain the very approximate 
expression, 






V 

18, 



-^( 4 V + 7 *iei + 4^ + 4^ + 7*2*2 + 4^) 
~ 360 ( 4 ^ 2 + 7 % e 3 + 4 e 3 2 ) - -^ (* 3 + e 3 ) - -^^j- 



y> (144.) 



945 



945 



5" ("l 2 + TB "l e l + 6 1 2 + ^ + 16 ^2 *2 + *2 2 ) 

V'3 T 16 "3 e 3 "I- e 3 y 40320 725760 



J 



which is accordingly the sum of the terms of the fourth and sixth dimensions in the 
development of the rigorous expression (120.), and gives, by our general method, 
correspondingly approximate expressions for the integrals of disturbed motion, under 
the forms 



8 Sj 8 S 2 SS; 

H + h + 8,, 



'3 

8«,' 



TSn 



SSj 8Sg , « " 3 

8>j 2 "T 8» "•" 8, 



8S, 



8>J 2 



"3- 8,3+ 8,* + 



8S_ 3 
la' 
8S_ 8 
8>l 3 ' 



(145.) 



and 



8S, 
Pi = -17, 



P2 



8<? 



8S, 



8S, 8S 2 



8S 3 

8S_ 3 
8e 2 

8S„ 



^3— 8e Q 8e a ~"8< 



(146.) 



28. To illustrate by the same example the theory of gradually varying elements, 
let us establish the following definitions, for the present disturbed motion, 

h ~ W 1J X 2 — ™2> h ~ m 3 + g f > J 

and let us call these six quantities x, x » 2 « 3 X x X 2 X 3 the varying elements of that motion, 
by analogy to the six constant quantities e l e 2 e 3 p l p 2 p z , which may, for the un- 
disturbed motion, be represented in a similar way, namely, by (127.) and (128.), 
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e x = ri x — ^ t, e 2 = n 2 - 


■ ™2 *> e 3 


P\ ~ W 13 P2 — ™2> 


P3 






(148.) 



We shall then have rigorously, for the six disturbed variables n x ?? 2 r\ 3 ^i sr 2 * 3 , ex- 
pressions of the same forms as in the integrals (127.) and (128.) of undisturbed 
motion, but with variable instead of constant elements, namely, the following : 



*1 


= *1 + 


*1 


t, 


*2 


= 


%2 


+ 


h 


*, 


% 


w l 


= x 1? 






m 2 


— 


^•23 








sr 3 



= *? + h * — -g 8 12 > 
= X 3 -gt; 



(149.) 



and the rigorous determination of the six varying elements z x » 2 x 3 \ x k 2 K 3 , as func- 
tions of the time and of their own initial values e x e 2 e z p x p 2 p 3 , depends on the in- 
tegration of the 6 following equations, in ordinary differentials of the first order, of 
the forms (105.) : 



= fi; = v * (*i + x i *)> 



dt 

dx, 



SH 



and 



dt 
dt 



d\ x 
~di 



2 " "3 „2 



Sx 2 



y? t («2 + X 2 #), 



= ^§ = > 2 '(*3 + M-j^ 2 ), 



(150.) 



dK 



9 



8H, 



S« 



= — f* 2 (*i + h t), 
= —p 2 (* 2 + x 2 *), 



II 2 being here the expression 

H 2 = ^ {(*! + V) 2 + (* 2 + M) 2 } + |K*3 + M - l£' 2 ) 2 > 



(151.) 



(152.) 



which is obtained from (125.) by substituting for the disturbed coordinates n x n 2 % 
their values (149.), as functions of the varying elements and of the time. It is not 
difficult to integrate rigorously this system of equations (150.) and (151.) ; and we 
shall soon have occasion to state their complete and accurate integrals : but we shall 
continue for a while to treat these rigorous integrals as unknown, that we may take 
this opportunity to exemplify our general method of indefinite approximation, for all 
such dynamical questions, founded on the properties of the functions of elements C 
and E. Of these two functions either may be employed, and we shall use here the 
function C. 

29. This function, by (109.) and (152.), may rigorously be expressed as follows : 
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+£/*{(v-^y-v}^j 

and has therefore the following for a first approximate value, obtained by treating 
the elements x x x 2 x s X t X 2 X 3 as constant and equal to their initial values e x e 2 e 3 p x p 2 p^, 

C = ~ 1 {/** W + ^) + " 2 ^ } + i {l* 2 (Pi 2 + ft 2 ) + * 2 ft 2 } 1 



--^gft + ^'V' 



In like manner we have, as first approximations, of the kind expressed by the ge- 
neral formula (Z 1 ), the following results deduced from the equations (151.), 



X 2 = p 2 - y? \e 2 t + \ P 2 t 2 ), 

h^Pz-^i^t + ^PsP-^gt 3 ), _ 
and therefore, as approximations of the same kind, 



(155.) 



e i = - •« Pi * 



e 2 — '■ ~ o ft ^ ~ 



/St > 



Pst + ^gt 2 - 



■Pa 



(156.) 



: 3-~-2«'T-p« — v 2jf • 

Substituting these values for the initial constants e x e 2 e 3 in the approximate value 
(154.) for the function of elements C, we obtain the following approximate expres- 
sion C x for that function, of the form supposed by our theory : 



Ci 



_ 1 f ih -p^ ± ^ -^ i. ( ^ 8 ~ ^' 



2 / 



{ 



r 



+ 



} 



- 4 { ( X l "" Pl>Pl + ( X 2 ~ ft) Pi + (*3 ~ ft) (ft ~ 4 # } > ( 157> ) 
+ 24 J ^ 2 (ft 2 + ft 2 ) + " 2 ft 2 J - 24 " 2 5" ft + 90 " 2 ^ 2 - 

The rigorous function C must satisfy, in the present question, by the principles of the 
eighteenth number, the partial differential equation, 

SC u? f/8C , V /8C • \ 2 ~) v s AC , 1 „\2 , 

and if it be put under the form (U 1 .), 

C = C x + C 2 , 
Cj being a first approximation, supposed to vanish with the time, then the correction 
C 2 must satisfy rigorously the condition 
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c=x'{ -v? +m + ^>y+^ + ^y+m^'-^y}*'] 



uv(w+^y+°^)><- 



f(l59.) 



In passing to a second approximation we may neglect the second definite integral, 
and may calculate the first by the help of the approximate equations (155.) ; which 
give, in this manner, 

C * = 'SI { (h ~ J°i) 2 + ( X 2 ~ ft) 2 + (h ~ ft) 2 } dt 

+ ifo* (^ - 4 s') ^ ~ ft) t2 d f 



= - ~ {(h - Pl f + (h - ft) 2 + <h ~ ft) 2 > 



(160.) 



+ 24 *f* 2 ft (*1 - ft) + f* 2 ft (*2 - ft) + " 2 ft (h - ft) ) 
~ 5 *S (h - ft) + 4 (^ft 2 + ^ft 2 + V >3 2 ) 

-Wo^SPz + ^S 2 - 

We might improve this second approximation in like manner, by calculating a new 
definite integral C 3 , with the help of the following more approximate forms for the 
relations between the varying elements X x X 2 K 3 and the initial constants, deduced by 
our general method : 

1 ~ 8^1 ~ tPi ~ ~ V?t V 1 "T 6 "T - 24 / "~ 2 V 1 "T" i« T 60 /> 

8C, 8C 2 _ A 2 -j? 2 / ^ f^\ fj> g / ^» ^\ J 

e 2 — - *p s ~ 8p 2 — ~~ ^ V 1 "T" 6 "^ 24 / "" 2 \ X "T" 12 "*" 60 /' 

„ - !£i ic, _ v-^3 (, , ££ , ^\ ^§ /, , v 3 j! , tf\ 16L) 

3 ~~ 8 ^a~ ^3 — ~ y 3 * V 1 "T 6 "•" 24/ "" 2 V 1 "•" 12 "•" 60/ 
"T 6 I 1 " 1 " 60 + 40/ ; 

in which we can only depend on the terms as far as the second order, but which 
acquire a correctness of the fourth order when cleared of the small divisors, and give 
then 



h=Pi-t* 2 t (ei + \ Pi*) + y ^ & ( e i + \P\*)> 

h = ft - V?t (e 2 + a" ft V + "5" ^ 4 * 3 ( e 2 + Tft V' 

^3 = ft-^K e 3+¥ft^-¥^^)+^ v4;3 G3+-4-ft ; -i^ #2 )' 



> . (162.) 
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(163.) 



But a little attention to the nature of this process shows that all the successive cor- 
rections to which it conducts can be only rational and integer and homogeneous 
functions, of the second dimension, of the quantities X x X 2 K 3 p l p 2 p 3 g, and that they 
may all be put under the following form, which is therefore the form of their sum, 
or of the whole sought function C ; 

C = p ~ 2 % (h - Pl ) 2 + b /tPl (K t - Pl ) + p* CfiPl 2 
+ ^~ 2 % (h ~ P2) 2 + KP2 (h-Pz) + I* ^ft 2 

+ "~ 2 «v (^3 ~ PsY + h *Pl (h - Pi) + » 2 C, p A 2 

+ f,g (h - Ps) + " 2 KgPz + » 2 ig 2 , 
the coefficients a^ a, &c. being functions of the small quantities p, v, and also of the 
time, of which it remains to discover the forms. Denoting therefore their differen- 
tials, taken with respect to the time, as follows, 

ddp = a'pdt, da f = a\dt, &c, (164.) 

and substituting the expression (163.) in the rigorous partial differential equation 
(158.), we are conducted to the six following equations in ordinary differentials of 
the first order : 
2«',= (2a„+^ 2 2 ; b',= (2a, + **t)(b, + t)>, d, = \ (6, + *) 2 ; 



/»,« (2«> + i»0 (/„- f * 2 ); K= (ft, +*)(/.- \fi); *',=!(/, 



' = 4 ( f — i 



* 2 ) 2 ;. 



(165.) 



along with the 6 following conditions, to determine the 6 arbitrary constants intro- 
duced by integration, 



c 



£. it • ~t- 

24 ; "0— — 245 h —90- 



(166.) 



a o 2 1 » "0 — 2 ' •/ ~" 6 

In this manner we find, without difficulty, observing that a^ b^ c^ may be formed 
from «„ ft, c, by changing v to p, 



ff. 



-i»a* 



v cotan vt, «« = — 5 fi 2 ^ ~ i j"- cotan /* t, 



&, = - ' + T tan 2"' 

c » = - 27 + ^ tan "2 ' 
/ = |^--i+-f cotan^, 

^v — 2 y 3 "" V 3 n 2 ' 



~ 6V - 2^ COtan ' *' 



*,= 



*+-tan^ 



2 ' 



"~ a m 2 ■+" (* 3 tan 2 ' 



(167.) 



The form of the function C is therefore entirely known, and we have for this/unc- 
tion of elements the following rigorous expression, 
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(\~Pif + (h ~ PzY (h - Psf 



2 m tan jw. t 



2 y tan v t 



> (168.) 



- T ^ - ^) 2 + ( A 2 ~ ft) 2 + (*s ~ ft) 2 } 

- t {ft (Xj - ft) + ft (X 2 - ft) + ft (X 3 - ft) } 

+ — {ft (Xj - ft) + ft (X 2 - ft)} tan £$■ + — ft (X 3 - ft) tan ~ 
~ j (ft 2 + ft 2 + ft 2 ) + ~ (ft 2 + ft 2 ) tan & + -f ft2 tan ^ 
+ (y ~ "Jr +4cotanv*)g(X 3 -ft)+ (-|- - -f tan ^) gft 

+ W ~ ¥ ~~ 2~v cotan " # 2 > 

which may be variously transformed, and gives by our general method the following 
systems of rigorous integrals of the differential equations of varying elements, (150.), 
(151.): 



P - _ l£ - _ A *--Pi _ Si tan £-' 



e 2 = 



8_C 



VZZ2 _ A tan £? 

jK.sin(*« ft lctu 2 ' 



>. . (169.) 



and 



_ £> tan - 4- &■ (— _ 1\ 

lp 3 ~ vsmvt v 2 ~ v Vsinvif y /» 

* x = g-^ = - (X x - ft) (* + j cotanp*) +ft (- t + -i- tan ^), 

*2 = st;=- (*2-ft)(* + 7 cotan p + ft (~^ + 7 tan lD> 
*3 = 8l£ = - (*3 ~ ft) V + 7" cotan ' +ft (~ * + T tan ^) 
+ ^("2"-"^ + 7" C0tan "V ; 



>(170.) 



that is, 



X a = ft cos jMi ^ — e x \j> sin ^ tf, 
Xg = ft cos pi — e 2 p sin ^ £, 

X 3 = ft cos v t — e 3 c sin v t + g" (^ sin c £ ), { 



1 



071.) 



and 



»j = e x (cos j«( t + ju< £ sin p t) + ft ( — sin p t — t cos ,«- £), 

« 2 = e 2 (cos \h t + f« t sin p £) + ft (— sin ju» £ — £ cos ^ t), 

\ f* / 

« 3 = e 3 (cos v £ + v t sin » t) -\- ft (•— sin v t — £ cos v t) 



-s\r? — T smpt + Y)- 



(172.) 
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Accordingly, these rigorous expressions for the 6 varying elements, in the present 
dynamical question, agree with the results obtained by the ordinary methods of inte- 
gration from the 6 ordinary differential equations (150.) and (151.), and with those 
obtained by elimination from the equations (113.) (114.) (147.). 

Remarks on the foregoing Example. 

30. The example which has occupied us in the last six numbers is not altogether 
ideal, but is realised to some extent by the motion of a projectile in a void. For if 
we consider the earth as a sphere, of radius R, and suppose the accelerating force of 
gravity to vary inversely as the square of the distance r from its centre, and to be 

= g at the surface, this force will be represented generally by ^-5- ; and to adapt 

the differential equations (78.) to the motion of a projectile in a void, it will be suffi- 
cient to make 

U = g W (i - j9 ........... . ■ (173.) 

If we place the origin of rectangular coordinates at the earth's surface, and sup- 
pose the semiaxis of + 2 to be directed vertically upwards, we shall have 

r = V(R + *) 2 + * 2 + y\ ........... (174.) 

and 

v = - t . +*£-*&£& aw 

neglecting only those very small terms which have the square of the earth's radius 
for a divisor : neglecting therefore such terms, the force-function U in this question 
is of that form (110.) on which all the reasonings of the example have been founded ; 

the small constants p, v, being the real and imaginary quantities \/&, \/ -~ir*> 

respectively. We may therefore apply the results of the recent numbers to the 
motions of projectiles in a void, by substituting these values for the constants, and 
altering, where necessary, trigonometrical to exponential functions. But besides the 
theoretical facility and the little practical importance of researches respecting such 
projectiles, the results would only be accurate as far as the first negative power (in- 
clusive) of the earth's radius, because the expression (110.) for the force-function U 
is only accurate so far ; and therefore the rigorous and approximate investigations of 
the six preceding numbers, founded on that expression, are offered only as mathe- 
matical illustrations of a general method, extending to all problems of dynamics, at 
least to all those to which the law of living forces applies. 

Attracting Systems resumed : Differential Equations of internal or Relative Motion ; 

Integration, by the Principal Function. 

31. Returning now from this digression on the motion of a single point, to the 
more important study of an attracting or repelling system, let us resume the differen- 
tial equations (A.), which may be thus summed up : 
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dtlll = 2(«Mw — dw&jj); . (A 2 .) 

and in order to separate the absolute motion of the whole system in space from 
the motions of its points among themselves, let us choose the following marks of 

position : 

_ 2.m x 2 .my 2 ,m % 

X » ~ 2m > &» ~~ 2m > Z " 
and 



2m 



(176.) 



% = x i -w n ,fi i —y i — y n ,Z A = % i -z n ; (177.) 

that is, the 3 rectangular coordinates of the centre of gravity of the system, referred 
to an origin fixed in space, and the 3 n— 3 rectangular coordinates of the n— 1 masses 
m x m 2 . . . m n _ l3 referred to the nth mass m n , as an internal and moveable origin, but 
to axes parallel to the former. We then find, as in the former Essay, 

+ i2 / .m(r2 + V 2 + C 2 )-2^{(2 / .^r) 2 +(2 / .mV) 2 + (2 ; .m0 2 },} ' ^^ 
the sign of summation 2, referring to the first n — 1 masses only ; and therefore, 

T = zTm{ fcr) + \w) + (jj) } 
If then we put for abridgement, 



x > _ _L £1 _ p _ s i- m ? " 
' m 8£' ? 2m , 

1_ 8_T , 2 l .mtl , 



" / »? S m' ^ 



.2m 



* / — m 8 ?' 



1 ST _ ■?,.»? 



(180.) 



we shall have the expression 

H = i W + y\f + *'„ 2 ) 2 m + 1 2, . m (*? + y « + &',*) ] 

+ dr « 2 ' • w *'<) 2 + ( 2 * • w ^) 2 + (2, • » *',) 2 > - u, 



(B*.) 



} • (c 2 .) 



of which the variation is to be compared with the following form of (A 2 .), 
dthB. = (dx u haf u — dx'^x^ + dy^y' H — di/^y,, + dz^rf u — dz'^zj 2 m • 

+ 2 r m(d&af, - dx , l h% + drfy' l -dy'fy+d®z' l -d& l l %0, 
in order to form, by our general process, 6 n differential equations of motion of the 
first order, between the 6 n quantities x n y n z n x' u yj z' u I r\ £ x [ { y\ z\ and the time t. In 
thus taking the variation of H, we are to remember that the force- function U de- 
pends only on the 3 n — 3 internal coordinates g ^ £, being of the form 

s 2 
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XJ = m n (m 1 f 1 + m 2 f 2 + ... + m n _ x /„_,) j 

+' *"i m 2 f 1>2 + m 1 m 3 / 1)3 + . . . + m n _ a m n _ 1 f n _ % „ _ „ J 



(D 2 -) 



in which /^ is a function of the distance of m £ from m n , andf i>k is a function of the 
distance of m 8 - from m h , such that their derived functions or first differential coeffi- 
cients, taken with respect to the distances, express the laws of mutual repulsion, being 
negative in the case of attraction ; and then we obtain, as we desired, two separate 
groups of equations, for the motion of the whole system of points in space, and for 
the motions of those points among themselves ; namely, first, the group 

dy^y'.d^dy'^O, > (181.) 

d % n = z' u dt, d z' u = 0, 
and secondly the group 

d % = \' r '/ + ~zr %.ma?,) dt, d *', 

d % = \^ l + ^^r mz ') dt > dz \= : ^-ji dt - _ 

The six differential equations of the first order, (181.), between x n y n z u x' ll y' u z' ls 
and t, contain the law of rectilinear and uniform motion of the centre of gravity of the 
system ; and the 6n — 6 equations of the same order, (182.), between the 6n — 6 
variables ly^x^ y\ z\ and the time, are forms for the differential equations of internal 
or relative motion. We might eliminate the 3 n— 3 auxiliary variables od t y\ z\ between 
these last equations, and so obtain the following other group of 3 n — 3 equations of 
the second order, involving only the relative coordinates and the time, 



) dt, 



m 0j; ' 



s — m 8g 

j> - 1 m 

4 ~~ m S*j 



l_HJ j_ 8U 



~ m.. / 8« ' 



> 



(183.) 



1*1 = 1™ + JL*. 

^> m. h " ' m I 



8, 

8U 
m 8? "•" m n ~> 8? 

but it is better for many purposes to retain them under the forms (182.), omitting, 
however, for simplicity, the lower accents of the auxiliary variables rf )y' \ z'„ because it 
is easy to prove that these auxiliary variables (180.) are the components of centrobaric 
velocity, and because, in investigating the properties of internal or relative motion, 
we are at liberty to suppose that the centre of gravity of the system is fixed in space, 
at the origin of xy z. We may also, for simplicity, omit the lower accent of 2„ un- 
derstanding that the summations are to fall only on the first n — 1 masses, and de- 
noting for greater distinctness the nth mass by a separate symbol M ; and then we 
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may comprise the differential equations of relative motion in the following simplified 
formula, 

dtlll = '2.m(d&a! — dx f U + drt!/ — dy , *ti + dZlsf-dxW3, . (E 2 .) 
in which 

H = |2.m(^ + y2 + ^2) + _i r {(2.m<) 2 + (2.my)2 + (2.m<) 2 } - U. (F 2 .) 

And the integrals of these equations of relative motion are contained (by our gene- 
ral method) in the formula 

SS = 2. m(a/&g — a'&a+y&u - 4'&j3 + *'&£— dly), . . (G«.) 

in which a |3 7 a' V c' denote the initial values of § n £ a?' y' 2', and S is the principal 

function of relative motion of the system ; that is, the former function S, simplified by 

the omission of the part which vanishes when the centre of gravity is fixed, and which 

gives in general the laws of motion of that centre, or the integrals of the equations (181 .). 

Second Example : Case of a Ternary or Multiple System with one Predominant Mass ; 
Equations of the undisturbed motions of the other masses about this, in their seve- 
ral Binary Systems ; Differentials of all their Elements, expressed by the coeffi- 
cients of one Disturbing Function. 

32. Let us now suppose that the n — 1 masses m are small in comparison with the 
nth mass M ; and let us separate the expression (F 2 .) for H into the two following parts, 

H x = 2.f (i + 5) (x' 2 +y' 2 + *'*)- M2. mf 



H 2 = 



M 



{x\ x' 2 + y\y\ + a\ s' 2 - M/ 1>2 ) + . . . > 



+ 



m,i m k 



M 



-(^+y,y*+^-M/ 44 ) + .... 



(H 2 .) 



of which the latter is small in comparison with the former, and may be neglected in 
a first approximation. Suppressing it accordingly, we are conducted to the following 
Qn — 6 differential equations of the 1st order, belonging to a simpler motion, which 
may be called the undisturbed : 

dt ~ 



m 8£ — in 8£' 



dt ~ m 8y — V ^ UJ W ' 

dt— m 83/' — V 1 "*" M/# ' dt "~ m 8 q ~~ m 8 q ' 

d< — 01 8s' — V 1 + M/ z ' 






"J 
H, 



1 
8/ 



At 8{f — M 8£ 



(P.) 



These equations arrange themselves in n — 1 groups, corresponding to the n — 1 
binary systems (m, M) ; and it is easy to integrate the equations of each gi'oup sepa- 
rately. We may suppose, then, these integrals found, under the forms, 
k = tf» (*, g, % I, xf, y>, z% v = % (4) (*, |, „, 5, *', y, *'),] 

X = tf» (t, g, n, £ at, y', z>), r = tf*) (t, g, „, & a/, 3/', z'), \ . . (K 2 .) 

i» = * (3) & i, 9, & ^ y, «'), » = ^ 6) (*, s, *, & ^ 3/', «o J 
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the six quantities kX^vtu being constant for the undisturbed motion of any one 
binary system ; and therefore the six functions yp), y^\ %( 3 ), %W, %( s ), %^>, or «, X, y*, 
v, <r, a, being such as to satisfy identically the following equation, 

Sx . 8x8 Hi 



° — m S* + S£ 8*' 



^8Hi 8xS_Hi Sx 8H t 8x8H t 
+ 8^ 8y ~"8y 8ij + 



_S_x SH, 
8 a' 



(LI) 



8*» ""8f"'~ r 8i) 8y 8y 8ij T S? Sz' 8a' 8?' 

with five other equations analogous, for the five other elements X, p, v, r, a>, in any one 
binary system (to, M). 

33. Returning now to the original multiple system, we may retain as definitions 
the equations (K 2 .), but then we can no longer consider the elements x. X. p. v. r. &>. of 

the binary system (m., M) as constant, because this system is now disturbed by the 
other masses m h ; however, the 6 n — 6 equations of disturbed relative motion, when 
put under the forms 

drf SJHi SJHg 

dt 
da' 



m 



m 



dt 

dr\ 
It 



_ 8 H x 8 H 2 
~" 8^ "•" 8^ ' 

8 H. 8 H 2 

"' sy + &y > 



m 



dt 

dy' 



SH, 



m-rr — 



H 

Si ' 



> 



*1 
dt 



SH, . 8H, 



m H7 = TJ + 






8JHj 

8? 



8H a 



(M 2 .) 



8 a" "* d * ~~ " 8 ? ' " 8 ? ' 

and combined with the identical equations of the kind (B 2 .), give the following simple 
expression for the differential of the element #, in its disturbed and variable state, 



m 



dx 
"dt 



Sx 8 H 3 Sx SHg Sx^ hU % 



S£ 8*' 8^' 8£ ~ 8>j 8y 8y Sij ~ 8? 8a' 8^ 8? 5 



k8H 2 Sx 8H 2 



Sx SH, 



(N 2 .) 



together with analogous expressions for the differentials of the other elements. And 
if we express \n Z,a/ y' z', and therefore H 2 itself, as depending on the time and on 
these varying elements, we may transform the 6 n — 6 differential equations of the 
1st order, (M 2 .), between % n Z, x' y* z' t, into the same number of equations of the 
same order between the varying elements and the time ; which will be of the forms 



TO 



dx 
It 



8H, 



8H, 



SH, 



8H, 



8H, 1 



{*, M Tx"+ {x ^ } "87 + {*> " } TV + {*» T > T7 + <*» "> ~sV 



dx 8H 2 



SH, 



8H, 



SH, 



8H, 



d£ 



dp 



TO 



V- 



Sx 
8H 



+ {*»**> -nr + <*> »> T7- + ^ T >T7 ' + <*» ^XT' 



8H, 



8H, 



SH, 



^7 = it*,*} Yf + {[*,-k}-i^+ {p, >}~ff+ It*,*} -jf+ ^*>T 



8H, 



m 



m 



dv 
d* 

dr 



8H, 



SH, 



SH, 



SH, 



{', •> Tf+ l'> x > Tx + <'» ^> 17T+ <"> T > TT + <'» •> 7 



H 



h° 2 o 



Sx 

SH 



8H 2 r _ „ 1 S_H 2 , r _ 8H, f _ JH, 



rfT = (*> *) TT + <«■» X >T7T + ^f > V+ {r ' v} TT+ {r ' *>TT* 



d 



8H 



TO 



8H 



"d? = t»* «} "FT + {*>> ^} TF" + ^»/*> 



8 Ho ,8 Ho . . .8 H, 
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if we put, for abridgement, 



r i—iifi? 1 Sx SX Sx SA Sx8a 8_x Sx 



8 x 8a. 

8~z'8~?' 



(P 2 -) 



and form the other symbols {%,(*}, {k,x}, &c, from this, by interchanging the letters. 
It is evident that these symbols have the properties, 

{X, «} = - {«, k}, {»,*} = ; . . . . (184.) 

and it results from the principles of the 15th number, that these combinations {z, k}. 
&c, when expressed as functions of the elements, do not contain the time explicitly. 
There are in general, by (184.), only 15 such distinct combinations for each of the 
n — 1 binary systems ; but there would thus be, in all, 15 n — 15, if they admitted 
of no further reductions : however, it results from the principles of the 16th number, 
that 12 n — 12 of these combinations may be made to vanish by a suitable choice of 
the elements. The following is another way of effecting as great a simplification, at 
least for that extensive class of cases in which the undisturbed distance between the 
two points of each binary system (m, M) admits of a minimum value. 

Simplification of' the Differential Expressions by a suitable choice of the Elements, 

34. When the undisturbed distance r of m from M admits of such a minimum q, 
corresponding to a time r, and satisfying at that time the conditions 

r' = 0, r">0, , (185.) 

then the integrals of the group (I 2 .), or the known rules of the undisturbed motion of 
m about M, may be presented in the following manner : 

«= ^{(iy-^.') 2 + (^'~^y) 2 + (^'-^0 2 }; 

k = a — §y' -\-qaf; 






=<-x 



x/ 



M 



dr 



M + m V dr* 



. dr 



> ■ (Q 2 .) 



m ■=. v -]- sin 



X /{ 2f c + £M/(,)-(l + g)^}' 

/ M + m 



x.%1 — 



' ^2A 






dr 

Vdr* 



dr 



y/{s* +«*/■(»■) -(1+5)5}' 



the minimum distance q being a function of the two elements %, p, whicl must satisfy 
the conditions 

af» + 2M/( ? )-(l+5)^ = 0,M/(y) + (l+5)7>0i.- ( 186 -) 
and sin" * s, tan" 1 t, being used (according to Sir John Herschel's notation) to ex- 
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press, not the cosecant and cotangent, but the inverse functions corresponding to sine 
and cosine, or the arcs which are more commonly called arc (sin = s), arc (tan = t), 

d V 

It must also be observed that the factor ~7j=%, which we have introduced under the 

signs of integration, is not superfluous, but is designed to be taken as equal to posi- 
tive or negative unity, according as d r is positive or negative ; that is, according as 
r is increasing or diminishing, so as to make the element under each integral sign 
constantly positive. In general, it appears to be a useful rule, though not always 
followed by analysts, to employ the real radical symbol ^R only for positive quan- 

T 

tities, unless the negative sign be expressly prefixed ; and then — p= will denote posi- 
tive or negative unity, according as r is positive or negative. The arc given by its 
sine, in the expression of the element &>, is supposed to be so chosen as to increase 
continually with the time. 

35. After these remarks on the notation, let us apply the formula (P 2 .) to calculate 
the values of the 15 combinations such as {x, X}, of the 6 constants or elements (Q 2 .). 

Since 

r= s^ + ^ + Z*), (187.) 

it is easy to perceive that the six combinations of the 4 first elements are as follows : 
U,X} = 0, U,p} = 0, {«,»} = 0, {\,p} = 0, {X,i>} = 1, {[*,»} = 0. . (188.) 
To form the 4 combinations of these 4 first elements with r, we may observe, that 
this 5th element r, as expressed in (Q 2 .), involves explicitly (besides the time) the 
distance r, and the two elements st, p ; but the combinations already determined 
show that these two elements may be treated as constant in forming the four combi- 
nations now sought ; we need only attend, therefore, to the variation of r, and if we 
interpret by the rule (P 2 .) the symbols {*, r} {X,r} {p,r} {v,r}, and attend to the 

equations (I 2 .), we see that 

(lv 

{*, r} = 0, {X, r} = 0, {p, r} = - Tp {v,r} = 0, (189.) 

dv 

Y t being the total differential coefficient of r in the undisturbed motion, as determined 

by the equations (I 2 .) ; and, therefore, that 

{*,r} = 0, {X,r}=0, {e,r} = 0, ■ . . (190.) 

and 

8 t dr , dt dr . 
{[*,r}=- rrI} =+^ I - t = l: (191.) 

observing that in differentiating the expressions of the elements (Q 2 .), we may treat 
those elements as constant, if we change the differentials of In £ oo' y' z' to their undis- 
turbed values. It remains to calculate the 5 combinations of these 5 elements with 
the last element m ; which is given by (Q 2 .) as a function of the distance r, the coor- 
dinate Z,, and the 4 elements x,, X,(a,v; so that we may employ this formula, 

te»} = 17 {e,r} + jf {e, Z,} + ^ {e, *} + jj- {e, X} + -^ {e,p} + -^ {e, v), (192.) 
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in which, if e be any of the first five elements, or the distance r, 

1 / tie tie 8<?\ tie 

{e,r} = --{ZM+*jjf + Z,$j)>{e>Z > }==-Yj,{e,z}==0, . (193.) 

and 

If""W/ > 17 — ~drT$> 17 - 1 ' ( 194 v 

the formula (192.) may therefore be thus written: 



fe.)-( 'iJh!Jy + tfe) ..\(£)-> 



(195.) 



W CW CO 

+ fa '> + TX fa x > + V fa M 
We easily find, by this formula, that 

Sri. /ul — fl • /*• /•.■I — 



{*»*} = — 1; {X,«} = 0; {f*,»} = 0; {r,*} ss-jjyj; . . . (196.) 



and 

, 8v Sco 8« 

{"^} = -85F?-u = ° (197.) 

The formula (195.) extends to the combination {r, «} also ; but in calculating this 
last combination we are to remember that r is given by (Q 2 .) as a function of », p, r, 
such that 

8r __ dt 

tir~~dr ; (198.) 

and thus we see, with the help of the combinations (196.) already determined, that 

. , St Sco 8 /*" . 8 pr 

{rja , } = ___ 8 -- = :_y ? ® r dr + r ^J q Cl r dr, (199.) 

if we represent for abridgement by @ r and Cl r the coefficients of d r under the integral 
signs in (Q 2 .), namely, 



' . . (200.) 



O — x . / M ± m dr fo I o TV/T J?/ N M + TO x 3 •) 



-* 



These coefficients are evidently connected by the relation 

8© 8/2 

if + T7r :=0 ' • • • • (201.) 



which gives 

J 
8 



^X^rdr+^/J^dr^O, (202<) 

r t being any quantity which does not vary with the elements k and p ; we might 
therefore at once conclude by (199.) that the combination {r, a} vanishes, if a diffi- 
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culty were not occasioned by the necessity of varying the lower limit q, which de- 
pends on those two elements, and by the circumstance that at this lower limit the 
coefficients Q r CL r become infinite. However, the relation (202.) shows that we may 
express this combination {r, a>} as follows : 

. <«"»•> = -kfi' & rdr + ^/J' Q r dn . . . ■ (203.) 

r t being an auxiliary and arbitrary quantity, which cannot really affect the result, 
but may be made to facilitate the calculation ; or in other words, we may assign to 
the distance r any arbitrary value, not varying for infinitesimal variations of *, /*, 
which may assist in calculating the value of the expression (199.). We may there- 
fore suppose that the increase of distance r — g'm small, and corresponds to a small 
positive interval of time t — r, during which the distance r and its differential coeffi- 
cient r' are constantly increasing ; and then after the first moment r, the quantity 

®r = V (204.) 

will be constantly finite, positive, and decreasing, during the same interval, so that 
its integral must be greater than if it had constantly its final value ; that is, 

t-r-J t g r ® r dr>{r-q)<d r (205.) 

Hence, although @ r tends to infinity, yet (r — q) & r tends to zero, when by dimi- 
nishing the interval we make r tend to q ; and therefore the following difference 

J, a >dr-- w -yJ 9 ® r dr = - w -J q {-^- Y )® r dr, . . (206.) 

will also tend to 0, and so will also its partial differential coefficient of the first order, 
taken with respect to /*. We find therefore the following formula for {r, a), (re- 
membering that this combination has been shown to be independent of r,) 

{r,»}=ri q {lfjQ r dr + ^^f;e r dr}; . . . (207.) 

the sign r _ implying that the limit is to be taken to which the expression tends 
when r tends to q. In this last formula, as in (199.), the integral / © r dr may be 

considered as a known function of r, q, k, p, or simply of r, q, %, if \h be eliminated 
by the first, condition (186.) ;, and since it vanishes independently of %, when r = q, it 
may be thus denoted : 

J^ r e r dr = <p(r,q,x)-<p(q,q,x), , . (208.) 

the form of the function <p depending on the law of attraction or repulsion. This 
integral therefore, when considered as depending on % and p, by depending on » 
and q, need not be varied with respect to *, in calculating {r, &>} by (207-), because 
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its partial differential coefficient ij^j O r dr), obtained by treating q as constant, 
vanishes at the limit r = q ; nor need it be varied with respect to q, because, by (186.), 

»g M+m * gg__ , ,„nqv 

Sue" 1 " M 2 a V — ^w.) 

It may therefore be treated as constant, and we find at last 

{«■,»}=<>, (210.) 

the two terms (199.) or (203.) both tending to infinity when r tends to q, but always 
destroying each other. 

36. Collecting now our results, and presenting for greater clearness each combi- 
nation under the two forms in which it occurs when the order of the elements is 
changed, we have, for each binary system, the following thirty expressions : 

{*, X} = 0, {*, /*} = 0, {*, »} = 0, {*, r} = 0, {*, «} = - 1, -J 
{X, z} = 0, {X, p) = 0, {X, v} = I, {X, r} = 0, {X, a) = 0, 
{j»,«} = 0, {p,X} — 0, {[*,»} = 0, {p,r} '= 1, {(b,a} = 0, 
{v, z) = 0, {», X} = - 1, .{v, p} = 0, {», t} = 0, {», *} = 0, 
{r, z} = 0, {r, A} = 0, {r, /»} = ~ 1, {r, «} = 0, {r, *} = 0, 
{*>, z} = 1, {«, X} = 0, {&>, /*} = 0, {&>, c} = 0, {&>, r} = ; 
so that the three combinations 

{(*, T} {«, «} {X, «} 

are each equal to positive unity ; the three inverse combinations 

{r, p} {«, a} {v, X} 

are each equal to negative unity ; and all the others vanish. The six differential 
equations of the first order, for the 6 varying elements of any one binary system 
(m, M), are therefore, by (O 2 .), 



(R».) 



d a. 

m lTt : = 

dco 
m-jT = 
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m -Tt = 
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6 p, 



> 



m 



(S 2 -) 



and, if we still omit the variation of t, they may all be summed up in this form for 
the variation of H 2 , 

m 2 = 2. w»(f6'&r -*>&/* + «/ &* — *'&« + X'&*-*'&X), . . (T 2 .) 

which single formula enables us to derive all the 6 n— 6 differential equations of the 
first order, for all the varying elements of all the binary systems, from the variation 
or from the partial differential coefficients of a single quantity H 2 , expressed as a 
function of those elements. 

t2 
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If we choose to introduce into the expression (T 2 .), for h H 2 , the variation of the 
time t, we have only to change hr to S r — S t, because, by (Q 2 .), & t enters only so 
accompanied ; that is, t enters only under the form t — r,-, in the expressions of 
f» *>i Ki ^i V\ z 'i as functions of the time and of the elements ; we have, therefore, 

a=-2^=-S.*V; (211.) 

and since, by (H 2 .), (Q 2 .), 

U 1 = '2.m^, (212.) 

we find finally, 

£S - _ Hi an ^ 

dt — It . . {V .) 

This remarkable form for the differential of H l3 considered as a varying element, 
is general for all problems of dynamics. It may be deduced by the general method 
from the formulae of the 13th and 14th numbers, which give 

lUi = S 2 (IMi iJLi _ lEi *Jh\ -l _i_ lii v (U*i *J** - LMi *J**\ 1 
dt 8x, V 8i) 8sr Sot 8jj / " t " " ' "*" 8x 6re \ 8>j § OT 8ot §,, 7 I 

' **! 8* "*" 8x 2 S# -I - - - - -t- g^^ 8* 8* ' J 

z 1 x 2 ..z 6n being any 6 w elements of a system expressed as functions of the time and 
of the quantities n *t ; or more concisely by this special consideration, that H^ + H 2 is 
constant in the disturbed motion, and that in taking the first total differential coeffi- 
cient of H 2 with respect to the time, the elements may by (F 1 .) be treated as constant. 
It is also a remarkable corollary of the general principles just referred to, but one not 

difficult to verify, that the first partial differential coefficient ^- s of any element x s , 

taken with respect to the time, may be expressed as a function of the elements alone, 
not involving the time explicitly. 

On the essential distinction between the Systems of Varying Elements considered in this 
Essay and those hitherto employed by mathematicians. 

37. When we shall have integrated the differential equations of varying elements 
(S 2 .), we can then calculate the varying relative coordinates £ n Z,, for any binary sy- 
stem (m, M), by the rules of undisturbed motion, as expressed by the equations (I 2 .), 
(Q 2 .), or by the following connected formulas : 



§ = r (cos 6 + — sin (0 — v) sin v\ 
V = r /'sin 6 — — sin (6 — v) cos v J, 
Z = -7 x/2X*-^ sin (6 - ») : 



(V 2 .) 
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in which the distance r is determined as a function of the time t and of the elements 
r, x, p, by the 5th equation (Q 2 .), and in which 



/M + m dr x_ 

6 = *>+f q ,r -TTTTTT M + m ?T » ^ 

Sj ^2 [b + 2 M/(r) jj- . js-j- 

g being still the minimum of r, when the orbit is treated as constant, and being- still 
connected with the elements x, p, by the first equation of condition (186.). In astro- 
nomical language, M is the sun, m a planet, § n £ are the heliocentric rectangular co- 
ordinates, r is the radius vector, d the longitude in the orbit, a the longitude of the 
perihelion, v of the node, 6 — a is the true anomaly, 6 — v the argument of latitude, 
\h the constant part of the half square of undisturbed heliocentric velocity, diminished 
in the ratio of the sun's mass (M) to the sum (M + m) of masses of sun and planet, 

x is the double of the areal velocity diminished in the same ratio, — is the versed sine 

of the inclination of the orbit, q the perihelion distance, and r the time of perihelion 
passage. The law of attraction or repulsion is here left undetermined ; for Newton's 
law, {Ij is the sun's mass divided by the axis major of the orbit taken negatively, and 
x is the square root of the semiparameter, multiplied by the sun's mass, and divided 
by the square root of the sum of the masses of sun and planet. But the varying 
ellipse or other orbit, which the foregoing formulae require, differs essentially (though 
little) from that hitherto employed by astronomers : because it gives correctly the 
heliocentric coordinates, but not the heliocentric components of velocity, without dif- 
ferentiating the elements in the calculation ; and therefore does not touch, but cuts, 
(though under a very small angle,) the actual heliocentric orbit, described under the 
influence of all the disturbing forces. 

38. For it results from the foregoing theory, that if we differentiate the expressions 
(V 2 .) for the heliocentric coordinates, without differentiating the elements, and then 
assign to those new varying elements their values as functions of the time, obtained 
from the equations (S 2 .), and deduce the centrobaric components of velocity by the 
formulae (I 2 .), or by the following : 

.... Mg , _ Mrf . _M?_ , 01 „, 

00 ~ M + m' V —M + m' Z ~"M. + m'' •••••■ V^-) 

then these centrobaric components will be the same functions of the time and of the 
new varying elements which might be otherwise deduced by elimination from the in- 
tegrals (Q 2 .), and will represent rigorously (by the extension given in the theory to 
those last-mentioned integrals) the components of velocity of the disturbed planet m, 
relatively to the centre of gravity of the whole solar system. We chose, as more 
suitable to the general course of our method, that these centrobaric components of 
velocity should be the auxiliary variables to be combined with the heliocentric co- 
ordinates, and to have their disturbed values rigorously expressed by the formulae 
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of undisturbed motion ; but in making this choice it became necessary to modify 
these latter formulae, and to determine a varying orbit essentially distinct in theory 
(though little differing in practice) from that conceived so beautifully by Lagrange. 
The orbit which he imagined was more simply connected with the heliocentric mo- 
tion of a single planet, since it gave, for such heliocentric motion, the velocity as well 
as the position ; the orbit which we have chosen is perhaps more closely combined 
with the conception of a multiple system, moving about its common centre of gravity, 
and influenced in every part by the actions of all the rest. Whichever orbit shall be 
hereafter adopted by astronomers, they will remember that both are equally fit to 
represent the celestial appearances, if the numeric elements of either set be suitably 
determined by observation, and the elements of the other set of orbits be deduced 
from these by calculation. Meantime mathematicians will judge, whether in sacri- 
ficing a part of the simplicity of that geometrical conception on which the theories of 
Lagrange and Poisson are founded, a simplicity of another kind has not been intro- 
duced, which was wanting in those admirable theories; by our having succeeded in 
expressing rigorously the differentials of all our own new varying elements through 
the coefficients of a single function : whereas it has seemed necessary hitherto to em- 
ploy one function for the Earth disturbed by Venus, and another function for Venus 
disturbed by the Earth. 

Integration of the Simplified Equations, which determine the new varying Elements. 

39. The simplified differential equations of varying elements, (S 2 .), are of the same 
form as the equations (A.), and may be integrated in a similar manner. If we put, 
for abridgement, 

and interpret similarly the symbols (^, a, X), &c, we can easily assign the variations 
of the following 8 combinations, (r, x, v) {p, u, X) (p, x, v) (r, u, X) (r, a, v) (f*, x, X) 
(r, x, X) (p, a, v) ; namely, 

h (t, x, v) = 2 . m (t h p — r 5 p + z & co — x h a + v & X — v 5 X ) — H 2 %t, ' 

^ (l«(,ft),X) = 2 . m (/a- ^ r — |0. S r + ft» 5 « — uy S » + X S v — X S r) — H 2 S#, 

5 (fju, x,v) = 2 . m (|o< & t — \b 5 r + x & a — z ci w -f- v & X — p & X ) — H 2 § t, 

6 (r,a>,X) = 2 . m (r^jO. — r &f* + ^o^o — a ^ K H" ^o^"o — X&v) — H 2 &£, 

5 (r, &), v) = 2 . m (r i \h — t & (& + a>o & *o ~ u ^ K + v & ^ — 'o ^o) ~ H 2 & £, 

6 {frXj'k) =2 . m ((Jb S r — \h h r + z 5 a> — x 5 a> -\- X 5 v — X I v) — H 2 lit, 

§ (r,x,\) = 2 . t«(tSj«i — r Sp + »5<y — « S«y + X Bv — Xhv) — H 2 §£, 

5 ([*,ai, p) = 2 . m ((!A 5 r — |a S r + a &# — <y & * -f- v 5 X — c S Xq) — H 2 e$ £, . 

« X p * r <y being the initial values of the varying elements zX p v r a. If, then, 
we consider, for example, the first of these 8 combinations (t, z, v), as a function of 



>(Y 2 .) 
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all the 3 n — 3 elements ^ « 4 X^, and of their initial values ,«/ 0j£ u 0>i X 0>i , involving also 
in general the time explicitly, we shall have the following forms for the 6 n — 6 
rigorous integrals of the 6 n— 6 equations (S 2 .) : 

m i x i= 8^-(^*/"); m i*o ti = - j^r.(« , »**0; > (Z 2 .) 

i 0,t 

and in like manner we can deduce forms for the same rigorous integrals, from any 
one of the eight combinations (Y 2 .). The determination of all the varying elements 
would therefore be fully accomplished, if we could find the complete expression for 
any one of these 8 combinations. 

40. A first approximate expression for any one of them can be found from the form 
under which we have supposed H 2 to be put, namely, as a function of the elements 
and of the time, which may be thus denoted : 

H 2 = H 2 (t, x 1} Xj, p u »!, r u a> v ... «„_!, X n _i, ib n _ ls v n _ 1} r n _ 1} » n-1 ) ; . (As.) 

by changing in this function the varying elements to their initial values, and em- 
ploying the following approximate integrals of the equations (S 2 .), 
. l r*iH a ,. i /"SH 2 

* = * + !- / H T^dt,Z=:* -±f n -£*dt, >...... (B3.) 

v ' m*s o ox ' u m%J o oco 

For if we denote, for example, the first of the 8 combinations (Y 2 .) by G, so that 

G = {r, «, ,}, (».) 

we shall have, as a first approximate value, 

G^VTH^ + ^ + '.^-H,}^ (m) 

and after thus expressing G x as a function of the time, and of the initial elements, we 
can eliminate the initial quantities of the forms r * » , and introduce in their stead 
the final quantities \u u X, so as to obtain an expression for G x of the kind supposed 
in (Z 2 .), namely, a function of the time t, the varying elements f* a X, and their initial 
values jM-o &> X . An approximate expression thus found may be corrected by a pro- 
cess of that kind, which has often been employed in this Essay for other similar pur- 
poses. For the function G, or the combination (r, », v), must satisfy rigorously, by 
(Y 2 .) (A 3 .), the following partial differential equation : 
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n 8G .iW, lSG > 1 SG 1 SG HG \ _ g . 

and each of the other analogous functions or combinations (Y 2 .) must satisfy an 
analogous equation : if then we change G to G t + G 2 , and neglect the squares and 
products of the coefficients of the small correction G 2 , G x being a first approximation 
such as that already found, we are conducted, as a second approximation, on prin- 
ciples already explained, to the following expression for this correction G 2 : 

^ /" f SG i , tt / 1 SG, l SG! l SG, \} , . ,„„. 

which may be continually and indefinitely improved by a repetition of the same pro- 
cess of correction. We may therefore, theoretically, consider the problem as solved ; 
but it must remain for future consideration, and perhaps for actual trial, to determine 
which of all these various processes of successive and indefinite approximation, de- 
duced in the present Essay and in the former, as corollaries of one general Method, 
and as consequences of one central Idea, is best adapted for numeric application, and 
for the mathematical study of phenomena. 



